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FOREWORD 


a sampied-oata  control  system  is  a linear  systemf  Its  cairpling 
device  is  the  only  element  not  found  in  conventional  servomechanisms . 

The  sampling  device  may  be  represented  by  an  impulse  modulator  which  can 
be  dsscribcd  simply  and  conveniently  in  either  the  time  domain  or  the 
frequency  domain.  By  using  the  impulso  modulator  equivalent,  one  can 
analyse  the  whole  sampled- data  system  either  in  the  time  domain  or  in  the 
frequency  domain  by  methods  similar  to  those  used  on  conventional  systems. 
This  analysis  has  already  been  described^"  in  principle,  and  the  formal 
techniques  of  ar.a3.ysis  have  been  laid  out.  The  great  similarity  between 
the  analysis  of  sampled-data  systems  and  that  of  conventional  systems 
permits  the  translation  of  many  techniques  which  have  been  used  on  con- 
ventional-system analysis  into  terms  which  apply  to  sampled-data  systems. 
The  function  of  this  report  is  to  describe  the  translated  techniques  in 
detail.  The  report  is  divided  into  three  parts i (l)  a description  of 
signal.'?  and  operations  in  sampled-data  system  components,  (2)  techniques 
for  system  studies  usinjj  component  characteristics,  and  (3)  use  of 
simplified  versions  of  the  analysis  techniques  described  in  design 
procedures. 

Signed 

Robert  ¥.  Sit tier 
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ABSTRACT 


The  main  advantage  from  realizing  the  similarity  between 
continuous-data  systems  and  sampled-data  systems  cesses  in  design  problems. 
System  design  is  the  art  of  compromising  between  what  is  desired  in  a 
certain  system  and  what  is  physically  realizable.  Whereas  a specific 
analysis  problem  has  a unique  solution,  very  often  no  specific  design 
problem  exists  for  a given  situation  and  even  if  it  did,  there  would  exist 
many  solutions  rather  than  one.  Accordingly,  the  technique  of  system 
design  is  usually  a matter  of  cut- and- try  wherein  certain  designs  are 
visualized,  analyzed,  and  modified.  To  be  workable  this  cut-and-try 
process  must  have  short  steps  based  on  simple  approximate  analysis 
procedures. 

In  this  report  analysis  is  discussed  before  design.  First, 

the  components  in  s sampled-data  system  a^*e  described  both  in  the  time 

and  frequency  domains.  Second,  techniques  for  flow  graph  (block  diagram) 

manipulations  are  illustrated  for  reducing  a complicated  block  diagram 

to  a simple  cascaded  equivalent.  Third,  responses  of  a cascade  of 

sampled-data  system  components  are  characterized  by  pole-zero  locations 

and  partial -fraction  expansions  in  the  s-plane,  the  K.  x -plane  and  the 
_ _aY/2 

C ' -plane.  The  use  of  error  ccefficients  and  root-locus  techniques 
are  illustrated.  Finally,  design  techniques  are  discussed  after  the 
analytical  techniques  are  presented. 


iii 
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1.  DESCRIPTION  OF  SIGNALS  AND  OPERATIONS  IN  SAMPLED-DATA  SYSTEM  fYTWPQWPWTQ 
The  key  to  the  system  of  analysis  used  here  is  to  compare  sanpling 
with  impulse  modulation.  Having  made  this  comparison,  one  can  represent 
all  elements  in  a sampled-data  system  by  simple  equivalent  blocks  and  can 
manipulate  these  blocks  in  the  manner  commonly  used  in  the  study  of 
continuous-data  systems « The  impulse  modulator  idea  has  been  presented 
in  earlier  papers,  "^^ut  the  function  of  this  report  is  to  show  how  it 
is  used  in  system  design.  The  following  paragraph  presents  a short  summary 
of  the  characteristics  of  all  the  fundamental  types  of  linear  components 
which  appear  in  sampled-data  control  systems.  The  remainder  of  Section  1.1 
shews  in  a detailed  fashion  ho*  the  response  charateriBtics  of  these  components 
are  described  in  both  the  frequency  domain  and  the  time  domain. 

1.1  The  Elementary  Components  in  Saapled-Data  Systems 

The  linear  components  in  a sampled-data  system  are  of  four 
kinds:  (l»  the  impulse  modulator,  which  .has  continuous  input  and  sanpled 

output,  (2)  those  conventional  elements  which  have  continuous  inputs  and 
continuous  outputs,  (3)  those  components  which  have  sampled  inputs  and 
saiipled  outputs,  and  (U)  those  components  which  have  sampled  inputs  and 
continuous  outputs. 

1.11  The  Impulse  Modulator 

In  sampled-data  systems  there  are  two  kinds  of  signals:  sampled 

and  continuous.  The  process  of  sampling  a continuous  signal  is  the  process 
of  selecting  ordinates  at  regularly-spaced  instants.  The  sampled  signal 

is  in  reality  a set  of  selected  ordinates,  but  it  could  be  conceived  as 

> • — ' : 

Linvill,  W.K.  and  Salzer,  J.M.,  Analysis  of  Control  Systems  Involving 
Digital  Computers  (To  be  published  ir"  l.R.E.  Proceedings)- 
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the  limiting  case  of  a continuous  wave  obtained  by  modulating  a carrier  ..f 
sharp;  tall  pulses  by  the  wava  to  be  sampled.  There  is  a one-to-one 
correspondence  between  the  successive  samples  of  the  signal  and  the  areas 
of  the  successive  impulses  in  the  modulated  wave.  Figure  1 shows  the 
characteristics  of  the  impulse  modulator.  Treating  the  sampled  ware  as 
the  limiting  case  of  a continuous  wave  allows  one  to  think  in  terms  of 
conventional  time  functions  and  their  Laplace  transforms. 


(c)  a continuous  exponential  wave 


(d)  the  corresponding  sequence  of 
impulses 


Fig.  1 Characteristics  of  the  impulse  modulator 
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A3  Fig.  i indicates,  the  carrier  savo  H«nn+oii  k»  •*»  ~ r^rir 

of  unit  impulses  separated  by  T seconds.  The  impulse  modulator's  output 
is  called  9^  (t)  and  is  related  to  the  input  9^(t)  by  Eq.  1. 

©^(t)  - 9i(t)  x i(t)  (l) 

The  carrier  wave  of  the  impulse  modulator  can  be  resolved  into 
its  various  frequency  components . When  this  is  done  it  is  represented 
by  an  infinite  series  as  in  Eq.  2. 


i(t) 


Jk  Cl  t,  where 


(2) 


— *•  = LN/  i.  • 

Since  the  carrier  of  the  impulse  modulator  has  an  infinite  number  of  exponential 
components  instead  of  the  two  which  correspond  to  a single  sinusoidal  carrier, 
the  impulse  modulator  output  has  many  sidebands  instead  of  only  two  as  in  the 
case  of  a sinusoidal  carrier.  From  this  fact  and  the  fact  that  superposition 
of  signals  applies  to  the  impulse  modulator,  the  relation  between  the 
impulse  modulator  input  and  its  output  in  the  frequency  domain  will  be  inferred. 

If  a single  exponential  component  of  frequency  s and  amplitude  A 
is  applied  to  the  input  of  the  impulse  modulator,  then  at  the  output  there 
appears  a family  of  signals  all  of  amplitude  A/T  and  of  frequencies  3,  B*$C~Lp 
a-j-TI,  3*j2A,  a-j2fl,  ...  . Thus,  if  9^(t)  - A £-8^,  the  output  0^*(t)  is 
defined  by  Eq.  3» 

#,*(0  - ?,y%(s*in  h.  o) 

n«  -s» 

55:  “ 1 

In  this  report  impulse -modulated  time  functions  will  be  distinguished 
from  continuous  ones  by  the  asteriskC*),  thus  S.*(t)  is  the  impulse 
modulated  9. (t). 
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ls a more  complicated  signal  is  applied  to  uhs  input  Ox  x*n6  impulse 


modulator.  the  complex  amplitudes  of  the  various  frequency  components  can 
be  represented  byK^(s).  Each  input  component  will  result  in  an  infinite 
set  of  sidebands.  From  superposition  Eq.  k follows. 


»et*Cs) 


•jn.fl). 


(U) 


Equation  U indicates  that  the  transform  of  the  impulse-modulated 
output  signal  is  periodic  in  the  frequency  domain,  or  that  ^ (s)  ■ *8^  (s+jm./'L). 
This  fact  is  consistent  with  the  fact  that  a periodic  time  function  has 
a sampled  frequency  function.  Conversely,  a sanpled  time  function  should 
have  a periodic  frequency  function,  which  it  does . 

Refer  again  to  Fig.  1.  The  continuous  tins  signal  is  an  exponentially 
decaying  step  which  starts  at  t * 0.  Because  the  carrier  wave  is  a set  ul 
narrow  pulses  which  is  even  about  the  origin,  the  carrier  pulse  at  t » 0 
occurs  at  just  the  time  when  the  discontinuity  in  8^(t)  occurs.  The 
half  of  the  carrier  impulse  occurring  for  t>0  is  multiplied  by  A. 

The  half  of  the  t • 0 impulse  occurring  for  t <1  0 is  multiplied  by  0. 

Hence  the  first  output  pulse  area  is  ~ . Figure  1 shows  that  the  impulse 
modulator  output  is  a collection  of  impulses  each  of  which  occurs  at  t«nT. 

The  transform  of  a unit  impulse  occurring  at  the  origin  is  1.  The  transform 
of  a unit  impulse  occurring  at  t-mT  is  £ . The  transform  of  a whole 

set  of  impulses  is  a sum  of  the  individual  transforms.  Thus  the  transform 
of  the  sequence  of  exponentially  decaying  impulses  has  the  value  given  by 
Eq.  5. 


oo 

\ , ,-amTjk-smT  A 

1 1 

/ A & £ - ? 
m“$* 

■ A 

■j^-Cs+aJT  " 2 

L _l 

-« 

See  Appendix  A for  a detailed  discussion  of  this  point. 
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For  the  sane  example  the  same  result  can  b®  obtained  frcs 
application  of  Eq.  4.  For  0,  (t)  an  exponentially  decaying  step*%(s)  is 

A 

Whence  __  . 

* ^-ga-  a Ar  i il 

- s+a+jn/i  A|i  ✓ -(s^a)T  ? * ' ' 


s+& 


*/e\  - 1 ^ 
Cs;  - * > 


All  impulse  modulator  inputs  having  Laplace  transforms  which  are  rational 
functions  in  o can  be  resolved  into  partial  fractions,  and  each  term  can  be 
subjected  to  the  summation  of  Eq..  6.  Consequently,  if  *4^(s)  is  a rational 
function  in  s,  <0^(s)  is  a rational  function  of  . In  summary,  sampling 

is  analogous  to  impulse  modulation.  The  input  wave  of  the  impulse  modulator 
ccn  be  related  to  the  output  in  the  time  domain  as  in  Eq.  1.  The  impulse 
modulator  input  and  output  can  be  related  in  the  frequency  domain  ae  in 
Eq.  U-.  If  »0(,(s)  is  a rational  function  of  s,  .0^(3)  is  a rational  function 

of  esT. 

4Ht 

1.12  Continuous  Filters  and  Discrete  Filters 

Continuouc  filters  are  conventional  filters  which  have  continuous 
inputs  and  continuous  outputs.  They  are  familiar  to  engineers  and  are 
described  in  terms  of  their  impulse  response  h(t),  or  by  their  transfer 
function  K(s).  Linear  electric  networks,  linear  amplifiers,  and  conventional 
linear  servo  components  are  all  examples  of  continuous  filters. 

A discrete  linear  filter  is  a filter  which  accepts  a sampled  data 
input  and  transmits  a sampled  data  output  which  is  linearly  related  to  the 
input.  A discrete  filter  will  accept  input  data  only  at  the  sampling  instants 
and  will  transmit  an  output  only  at  the  sampling  instants.  Figure  2 shows 


** 


Precisely,  these  filters  should  be  called  continuous-signal  filters  and 
discrete-signal  filters.  The  modifiers  refer  to  the  nature  of  the  signal 
and  not  to  the  nature  of  the  filter  itself.  Since  all  filters  used  will 
be  of  the  lumped-parameter  type  and  since  frequent  interchange  between 
continuous  and  discrete  signals  will  be  made,  all  continuous-signal 
filters  will  be  referred  to  as  continuous  filters  and  ail  discrete-signal 
filters  will  be  referred  to  as  discrete  filters. 
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the  parallelism  between  a continuous  linear  filter  and  a discrete  linear  filter* 
A digital  computer  with  a linear  program  is  equivalent  to  a linear  discrete 
filter.  Another  example  is  a continuous  filter  with  an  impulse-modulated 
input  and  an  impulse  modulator  at  its  output.  This  example  is  so  prevalent 
that  it  is  described  in  detail  in  Sec.  2. 


»*(s}  -^(sWs)  »**(s)  - ^*(s)  K*(b) 


Fig.  2.  Comparison  between  a continuous  linear  filter  and  a discrete 
linear  filter. 

The  operation  of  a linear  discrete  filter  can  be  described  in  the 
time  domain  b 7 a difference  equation  such  as  Eq.  7*. 

9o*(n?)  - ho<Si*(nT)  ♦ \\*Rn-l)T~]  + h^Rn-a)  TJ  * — - h^*  [(n-k)Tj*  — , 

(7) 

where* 

S^CnT)  * an  inpulse  whose  area  is  the  output  sanple  at  t - nT. 

6^(nT)  - an  impulse  whose  area  is  the  input  sample  at  t ■ nT. 
h,  - the  constants  of  the  filter. 

j£ 
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A discrete  filter  can  be  described  equally  well  in  the  frequency 

-sT 

domain  by  transfer  functions  which  are  rational  functions  of  £ . The 

transfer  function  of  the  system  is  the  ratio  between  the  Laplace  transform 
of  the  output  and  the  Laplace  transform  of  the  input.  Since  if  the  Laplace 
transform  of  G /*(t)  is ^#(s) - then  the  Laplace  transform  of  6^*(t-T)  is 
^(s)€  , Eq.  7 can  be  expressed  in  the  transform  domain  as  in  Eq.  8. 

•©♦^*(8)  « hj©^*(s)  ♦ h^^T"*7  »©j|*(s)  ♦ hg£  2sT>©^(s)  ♦ ... 

-ksT  * (°) 

hk^  ^ (s)  ♦ ... 

**sT 

The  transfer  fvinction  in  this  case  becomes  a power  series  in  £ . For 

most  practical  cases  the  power  series  can  be  expressed  in  closed  form  and  the 

-sT 

transfer  function  may  be  expressed  as  a rational  function  of  £ having 

-sT 

only  a limited  nunfcer  of  poles  in  the  £ plane.  Techniques  for  expressing 
the  transfer  functions  in  various  forms  are  described  in  Sections  1.211  and 
1.212. 

Note  that  both  discrete  and  continuous  filters  can  be  described 
in  either  the  time  domain  or  the  frequency  domain.  While  a discrete  filter 
can  operate  or.  only  impulse-modulated  inputs,  a continuous  filter  can  operate 
on  both  continuous  inputs  and  on  impulse-modulated  inputs.  In  both  cases 
the  output  is  continuous.  To  find  the  Laplace  transform  of  the  output 
of  the  continuous  filter  when  the  input  is  impulse-modulated,  merely  multiply 
the  transform  of  the  input  by  the  transfer  function  of  the  filter. 

1 .13  Summary  of  Characteristics  of  Elementary  Components 
All  linear  sampled-data  systems  may  be  represented  by  systems 
composed  of  a few  elementary  building  blocks.  The  elements  are  the  impulse 
modulator,  continuous  linear  filters,  and  discrete  filters.  Both  types 
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of  filters  are  described  by  transfer  functions  in  the  frequency  domain  or 
by  unit  signal  responses  in  the  time  domain.  The  impulse  modulator  is 
linear,  but  is  like  an  amplitude  modulator  rather-  than  a filter,  and 
consequently  cannot  be  described  by  a conventional  transfer  function. 

It  can  be  easily  described  in  either  the  time  domain  or  in  the  frequency 
domain  in  a manner  analogous  to  the  way  any  amplitude  modulator  would  be 
described,  as  Eqs.  1 and  U shew.. 

1 . 2 Techniques  for  Har.d3.lng  Signals  in  Samplcd-Data  System  Components 
Thus  far  the  types  of  component?  encountered  in  sampled-data 
systems  have  been  described.  The  signals  are  either  sampled  or  continuous. 

This  section  is  devoted  to  a review  of  how  conventional  techniques  for 
describing  signals  in  continuous  systems  car.  be  applied  to  the  problem  of 
signals  in  sampled-data  systems. 

In  sampled-data  systems  three  classes  of  situations  arise  so  far 
as  signal  filters  are  concemeds 

(1)  Continuous  signals  are  applied  to  continuous  fillers, 

(2)  Impulse-modulated  signals  are  applied  to  discrete  filters, 

(3)  Impulse-modulated  signals  are  applied  to  continuous  filters* 
The  first  case  is  the  conventional  problem.  A continuous  filter  is  described 
precisely  by  a transfer  function  or  by  an  impulse  response.  The  filter's 
precise  response  from  any  input  may  be  obtained  by  multiplying  the  transfer 
function  by  the  input  transform  and  then  performing  an  inverse  transformation* 
Alternatively  it  can  be  obtained  by  convolving  the  impulse  response  wi+h  the 
input  time  function.  Though  these  precise  analysis  techniques  are  always 
applicable,  a system  designer  is  usually  more  interested  in  the  nature 

of  the  responses  to  each  of  several  typical  inputs  than  in  the  exact  form 
of  any  one  response;  consequently,  he  usually  describes  the  response 
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characteristics  in  less  precise  but  more  practically  meaningful  form  than  the 
exact  transfer  func  lion  provides  . Tvo  approxliDa  ts  wfiys  stg  often  Xiscci  for 
describing  response  characteristics  of  conventional  continuous  systems: 

fl)  in  terms  of  the  amplitude— puuo?  cziaracoei.  j.stiCo  of  the  transfer  function, 
(2)  in  terms  of  the  nature  of  the  impulse  response# 

The  form  of  the  amplitude-phase  curves  cf  the  transfer  function 


allows  one  to  make  certain  generalizations  about  the  nature  of  the  response. 
If  the  amplitude  is  approximately  constant  within  the  pass-band  and  if  the 
phase  is  linear,  then  the  response  will  be  satisfactory  for  all  signals 
whose  frequency  band-width  is  less  than  the  cut-off  frequency#  The  rise- tine 
of  the  response  to  a step  is  approximately  equal  to  one-half  the  reciprocal 
of  the  band-width  in  cycles  per  second,  and  the  delay  is  proportional  to  the 
phase  slope.  How  much  deviation  from  constant  amplitude  and  linear  phase 
characteristics  can  be  tolerated  cannot  be  clearly  stated  in  general,  but 
rule-of-thumb  allows  a system  transfer  function  to  have  a resonant  peak 
in  the  response  of  about  30  present  above  the  average  value  in  the  pass 
band  (Sp  - 1.3)  • For  more  precise  conditions  the  engineer  usually  makes 
exact  calculation  of  the  transient  responses. 

The  supeiposition  integral  generalises  the  system  impulse  response 
to  describe  the  system  response  to  any  input.  Though  the  precise  form  of 
the  impulse  response  is  not  important,  its  nature  is  most  important.  Of 
importance  is  the  duration  cf  the  impulse  response  and  the  degree  cf 
oscillation  in  it.  The  nature  of  the  impulse  response  is  easily  correlated 
with  the  poles  and  residues  of  the  transfer  function.  Many  practical  systems 
have  only  one  pair  of  dominant  poles,  which  usually  are  complex  conjugates. 
The  amount  jind  duration  of  oscillation  in  the  impulse  response  is  usually 
determined  by  the  ratio  of  the  real  part  of  the  dominant  poles  to  their 

^ Brown,  O.S.,  and  Campbell,  D.P.,  Principles  of  Servomechanisms,  John  Wiley 
and  Sens,  Inc.,  New  York  19k 3 , 
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magnitude.  This  ratio  is  commonly  called  the  damping  ratio,  t A value 

of  £ in  the  neighborhood  of  0.?  is  usually  satisfactory.  Sines  most  transfer 

functions  of  conventional  systems  are  rational,  functions,  they  can  be 

described  in  terms  of  pole-zero  constellations.  For  continuous  SySwCmv  f 

simple  graphical  proc^Vires  exist  to  relate  the  transient  response  to  the 

pole-zero  pattern. These  procedures  are  not  described  here,  but  their 

analogues  are  described  in  connection  with  sampled-data  transient  response 

calculations.  The  preceding  few  sentences  have  summarized  the  techniques 

for  describing  signal  operations  in  continuous  systems.  The  analogues  to 

these  techniques  will  now  be  applied  in  discrete-signal  systems. 

1.21  Operations  on  Impulse-Modulated  Signals  by  Discrete  Filters 

As  was  stated  earlier,  a linear  system  having  sanpled  inputs  and 

sampled  outputs  usually  has  a transfer  function  which  is  a rational  function 
C-  sT 

Ui  . 


From  such  a transfer  function  it  is  always  possible  to  calculate 


precisely  the  output  corresponding  to  a given  input.  For  the  sanpled-data 
system  this  problem  is  actually  simpler  than  for  the  continuous  case,  but 
it  is  not  a problem  usually  encountered  by  engineers,  so  it  is  discussed 
in  detail. 

1.211  Precise  Description  of  Signals  in  Discrete  Filters 
In  the  same  sense  that  precise  description  of  response  of  continuous 
filtfers  can  be  obtained  from  their  transfer  functions  and  precise  description 

of  the  input  signal,  the  same  sort  of  description  can  be  made  of  responses 

“ST 

of  discrete  filters.  If  the  transfer  function  D(  6 ~0i)  is  a rational  function 

IT 


Ibid.,  p.  U9- 

Guillemin;.  3. A.,  Mathematics  of  Circuit  Analysis,  John  Wiley  and  Sons,  Inc. 
New  York. 

Truxal,  J.G.,  Servomechanism  Synthesis  through  Pole-Zero  Configurations, 
Sc.B.  Thesis,  M.I.T.,  1950. 


Report  R-222 


-11- 


of  £~Sl,  it 


may  be  written  as  a ratio  of  two  iKJlynomials  as  in  Eq.  9* 


s1?»  *£<*(s)  a ♦ a 

D(  t j - o_  ■ o lv 

# 


-31 


-2s? . 


— ♦ a 


' -msT 


L_  V 


7-7T 


(9) 


(a)  1 ♦ bx  £'ai  ♦ b^”"1*  — ♦ b*  6 " k 31 

From  this  transfer  function  (a)  may  be  expressed  in  terms  of  *9^  ”(3)  as 


in  Eq o 10 • 


»«0’(.)  ♦ bj o.o*(s)  rsT  * b^o*(S)  e‘2sT  ♦ — * \^0*c”ksT  - 

ao^#(s)  ♦ a-L*e*i*(s)C”sT  * “fi*<3)£"2ST* * ««*i*fB)£*,n8T. 


(10) 


Because  of  the  fact  that  each  term  in  Eq.  10  can  be  identified  immediately 

4*. 

as  the  transform  of  a time  function,  9 (nT)  may  be  expressed  in  terms  of 
O^CnT). 


«0*(nT)  ♦ b^jCn-ljTj  ♦ [Cn-2)^|  ♦ ♦ bk  eo*^n-k)Tj- 

».»i*  <"t>  * =i8i*  m * ♦ — * »n«i*  [»«>*} 


The  relation  in  Eq.  11  holds  for  all  values  of  n.  Given  a sequence  of  input 

samples,  ©^(nT),  and  a set  of  initial  conditions  of  8^  (nT)  one  can  use 

Eq.  11  to  find  8 (nT)  for  all  n-values  after  the  solution  is  started, 
o 

4*. 

For  a specific  example,  suppose  that  8^  (nT)  is  defined  for  all  values 

4^ 

from  n * -m  to  n — >»•  and  that  9q  (nT)  is  defined  for  all  values  from 
n - - k to  n ■ -1.  The  value  of  0 (nT)  may  be  evaluated  for  n - 0 by 
Eq.  11.  Using  the  computed  value,  the  value  of  0q  (0)  and  Eq.  11  again, 

4t 

one  can  find  9q  (T)„  Ey  repeating  the  process  one  can  develop  the  whole 
sequence  of  output  samples  which  result  from  the  given  input.  The  nature 
of  this  solution  is  identical  with  the  nature  of  the  solution  of  a differential 
equation  by  successive  integrations  to  obtain  an  integral  equation.  The 
number  of  initial  conditions  necessary  is  exactly  equal  to  the  number  of 
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transfer  function  poles  in  the  s -plana.  These  initial  conditions  correspond 
to  the  constsnts  of  integration  in  the  integral-equation  case.  Because 
finite  time  delays  »nd  not  di  fferentials  are  involved,  the  solution  of  Eq.  11 
is  much  simpler  than  that  of  an  integral  equation-  For  times  soon  after 
the  start  of  the  transient,  Eq.  11  is  very  convenient  to  use;  however,  it 
would  be  quite  tedious  to  obtain  the  output  for  a large  value  of  n from  this 
equation.  To  get  a picture  of  the  nature  of  the  transient  for  a longer  time, 
expand  the  transfer  function  of  Eq.  9 into  a partial  fraction  expansion  as 


in  Eq.  12. 


«rsT)  • ^ 


♦ — * 


e * 1 € -p5 


e‘81-e? 


may  be  related  to  t©^^(s)  as  in  Eq.  9« 
_ oA  ^*(b)  0(  n <X*(s) 

* <s>  “ ~ir ♦ --lir- — 

€ 


♦ * 


W._»©t*(s) 

' ft.  X 

f Si-B, 

O pk 


Equation  13  leads  one  to  break  up  the  output,  *0^  (s),  into  components, 
each  of  which  is  related  to  one  natural  mode  (one  pole)  of  the  system. 
The  system  of  Eq.  ll*  results. 


* (1) 


o 'c; 


- Pi 


£ - P, 


•0*  (s) 


€ - h 


A 

V£i  (8)  , 
—i — ?r 

1 *T  e 

•(n  m 

yr-^i  (s) , 


1’?Tt 

i 

°<V:  * 

i/^l  r a 1 

SZ2 ""  • 
i - s-  csl 
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Ths  equivalent,  set  of  time  functions  is  giver,  in  Eq.  1$. 


*(1) 

a * <Vi‘) 
”o  ' * 

. “ 1 
Pi 

v.  t.nT)  ♦ 

l 

1 

»T 

Jl) 

© 

o 

(n-l)T  , 

*(2) 

®o  &*) 

2 

w 

S^vnT)  ♦ 

1 

*T 

. *(?> 
o 

(n-l)T  , 

(15) 

xW 

e (nT) 

k 

0±*(nT ) - 

X 

*(k) 

© 

0 

(n-l)T  . 

jf. 

To  determine  initial  conditions  of  all  components  of  SQ  (n'f),  the  same 

* 

number  of  initial  values  of  G are  needed  as  when  Eq.  11  is  usad. 

o 

Tne  procedure  for  finding  the  initial  conditions  closely  parallels  the 
procedure  for  getting  initial  conditions  in  the  solution  of  conventional 
differential  equations.  If  all  values  of  p have  magnitude  great  er  than 
unity,  all  components  of  the  output  will  ultimately  die  out  with  no 
input.  If  any  value  of  p has  magnitude  equal  to  1,  its  component  will 
never  die  oat  regardless  of  what  the  input  does.  If  any  6 exists  having 
a magnitude  less  than  1,  the  response  will  grow  without  limit.  Other 
maniuplations  of  the  transfer  function  could  be  made  in  order  to  emphasize 
specific  aspects  of  the  response,  but  these  manipulations  will  not  be 
described  in  detail  at  this  point. 

This  pjcocedure  for  analyzing  a system  with  sampled  input  and 
sampled  output  is  analogous  to  the  operation  in  a continuous  system  of 
precisely  relating  an  output  time  function  to  an  input  time  function. 


* 


The  general  problem  of  relating  transforms  of  sampled  functions  to 
the  sampled  functions  will  be  treated  in  Appendix  B. 
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1,212  Approximate  Characterisation  of  operation  on  Signals 
in  Discrete  Filtoys 

The  nature  of  response  characteristics  of  c discrete  filter 
may  often  he  obtained  from  a study  of  the  amplitude  and  phase  of  the 
transfer  function  as  functions  of  the  complex  frequency  s.  To  see 
physically  the  connection  between  the  amplitude  and  phase  description 
of  this  SyS  t/Gm  and  that  of  a continuous  signal  system,  refer  to  Fig.  3. 


! 

carrier  1 

~T~~ ' 

oi(t) 

1 

JU- 

System  with  sampled 

e ft) 

0 

impulse 

Modulator 

CD 

r £ 

input  and  sampled 

•0  ^ Ideal 

^(s) 

i <(S) 

output  with  transfer 
function  F(  ^”3^) 

>©<  (s) 

Fig,  3.  A system  with  sampled  input  and  sampled  output  fed  by 
input,  sampler  and  with  an  ideal  low-pass  filter  at  the 
oxitput  to  act  as  a desampler. 

The  signal  Q^(s)  is  a conventional  continuous  signal  which  may  be  described 
in  the  frequency  domain  or  in  the  time  domain.  The  impulse  modulator 
samples  the  continuous  signals  and  thereby  adds  the  complementary  signals 
as  indicated  by  Eq.  U and  illustrated  by  Fig.  U.  As  the  signal  >9^  (»)  is 
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Figo  U.  Amplitude  spectra  of  a continuous  signal  and  of  the 
corresponding  impulse-modulated  signal. 


“*sT 

fed  into  the  system  having  a transfer  function  n(/_~  ),  each  component 

of  the  signal  is  handled  in  precisely  the  same  manner  as  each  other  component 
because  D(£  "')  is  periodic  in  frequency  of  period-*  i-  . One  can  either 

, -c? 

interpret  the  result  in  terms  of  the  variable  fc,  “ , in  which  case  the 

sanples  stand  out,  or  in  terms  of  the  variable  s,  in  which  case  the 

frequency  components  of  the  corresponding  unsampled  signal  stand  out.  If 

the  anplitude  and  phase  of  D(P,"S^)  are  plotted  as  a function  of  s for 

then  one  can  visualize  how  each  frequency  component  of  each  complementary 

component  of is  operated  upon.  Since  D(£.  S^)  is  a rational 
-sT 

function  of  ^ , it  is  very  easy  to  visualize  its  value  as  a quotient 

of  products  of  vectors  in  the  £ S^-plane.  Suppose  D(  8^)  has  a 
constellation  of  poles  and  zeros  as  shown  in  Fig.  5*  The  value  of  the 


amplitude  and  phase  functions  is  shown  in  Fig.  So . The  "pass  band* 

of  the  system  is  in  the  range  from  about  to  radians  per  second. 
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Fig.  5>.  The  plot  of  the  amplitude  said  phase  of  D(£””1)  correlated  with 
its  constellation  of  poles  and  seres  for 

D.  ? ( f °T  - 0-1) 

[CsT  - (1.1  - J0.3)]{£“S*  - (l.i-jo.3^ 

The  resonant  peak  occurs  at  the  frequency  for  which  the  point  on  the  unit 
circle  lies  nearest  the  pole  of  Q,  or  at  about  1 •ft/PO.  At  th???  points, 

also,  the  phase  has  the  steepest  slope. 

* 

To  relate*^  (s-*  to  (s),  merely  multiply  the  components  of 
(s)  whose  magnitudes  are  Indlc^w^d  in  Tig.  U by  the  complex  function  D 
plotted  in  Fig.  pb.  If  the  signal  ^ is  fed  into  an  ideal  lew-pass 
filter  which  rejects  all  the  cojuplenentary  signals  but  passes  the  pure 
signal  without-  ’delay  or  di3tort5.on,  the  overall  result  is  the  earae  as  if 


Report  R-222 


•17- 


the  unsampled  input  signal  s > were  fed  into  a c ontinuous  filter  having 
a transfer  function  equal  to  D in  tho  range  from  -.j-^/2  to  ‘ti-^/2  and 
equal  to  zero  outside  this  range. 

The  equivalent  to  tha  impulse  response  of  a continuous  system  is 
the  unit-sample  response  of  a sanpled-data  system.  The  unit-sample  response 
can  be  obtained  by  th«  technique  of  either  F>q.  11  or  Eq.  15 . From  a study 
of  the  poles  and  residues  of  D(£  ~‘i),  one  can  deduce  the  form  of  the 
unit-sample  response.  Just  as  in  the  continuous-signal  case,  if  the 
sampled-data  system  has  two  dominant  poles  (tvo  which  are  much  closer  to 
the  unit  circle  in  the  ( ‘-plane  than  all  other  poles),  one  can  characterize 

/J 


the  nature  of  the  response  in  terms  of  the 


(a) 


Fig.  6, 


Loci  of  constant  3 in  the  s-plane  and  their  corresponding  locations 
-sT 

in  the  £ -plane . 
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<P  - ,«T  i 

shows  loci  of  constant  ^ in  the  -plane.  The  value  of  has  the 

same  meaning  for  sampled  systems  as  it  does  for  continuous  systems.  It 

gives  a measure  ox  how  fast  an  oscillatory  transient  dies  out  relative 

to  its  period  of  oscillation.  For  values  of  £ much  smaller  than  0.7 

the  transients  are  highly  oscillatory#  for  values  of  J[  much  larger  than 

0.7  the  transients  are  highly  damped.  If  there  are  more  than  one  pair 

] of  dominant  polos.,  the  impulse  response  can  he  calculated  but  cannot  be 

predicted  readily  from  pole  positions  because  there  are  too  many  components. 

i 

I 

1 

j 


| 

i 

i 


I 
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2.  .RESPONSE  OF  CASCADED  STStmS 

In  the  preceding  section  methods  were  given  for  finding  the  response 
of  the  three  basic  types  of  sampled-data  system  components  - the  modulator* 
continuous -signal  filter  amu  discrete-signal  filter.  in  this  section  we 
treat  the  problem  of  finding  the  response  of  any  cascaded  sequence  of  such 
components . 

Evidently  the  calculation  of  overs? 1 response-  of  such  a system 
f-ould  be  done  stepwise  by  methods  already  given*  proceeding  from  component 
to  component.  However,  suen  a procedure  is  more  tedious  than  necessary  since 
any  combination  of  cascade  domponents  can  be  reduced  to  an  equivalent  simple 
form  with  the  same  response.  The  fundamental  equivalence  vhich  permits 
this  reduction  is  shown  in  Fig.  7*  The  reduction  consists  of  replacing 


(a) 


inputs 

only 


(b) 


Original  System 


e*(t) 

o- 


t 


e*(t) 

-o 


Equivalent  System 

Fig.  7*  Reduction  of  a cascade  system  to  a single  element, 
the  combination  of  a continuous  filter  followed  by  an  impulse  modulator  by 
an  equivalent  discrete  filter  when  the  inpit  signals  are  sampled. 

The  equivalence  is  derived  by  means  of  the  following  procedure. 
Apply  a sampled  input  6*(t).  The  transform  of  the  signal  applied  to  the 


impulse  modulator  is  then 
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^(s)  -icjcs)  . 5(s)  > (16) 

The  cutout  transform  is  then  calculated  from  this  result  by  the  summation 


•»  00 


procedure  of  Eq,  U* 

^(s)  = ~ f ♦ jn  ) - w 2 K 's  + •1n-t  + jnjTi). 

n a "oo  ^ * **^o  ^17) 

But  since  O^(t)  is  a purely  0 ample  a signal,  k©£(s)  is  periodic,  i.e., 

£ ♦ 

/Q^(s)  *»Ot(s  ♦ jnXl)  for  any  integer  n.  Therefore 


tQt 


7,1  =.o *r~i  , i t“ 

n 


£ 


K(s  ♦ jnjTi) 

Oo 


(18) 


so  that  the  transfer  function  of  an  equivalent  system  is 
K*(s)  - j ) K(s  ♦ in£\). 


VX7  / 


n - -»o 


1 

T. 


Since  = ^ _ K(s  * jn.fl)  is  a function  of  s which  is  multiplied  by 

"n~  ♦ - Oo 


“^(s;  to  give  Wp(s),  it  is  the  transfer  function  of  the  discrete  filter 

made  up  of  a continuous  filter  and  an  impulse  modulator.  Recall  that  the 

components  which  accept  sampled  signals  and  deliver  ScuupxSG  signals  have 

transfer  functions  which  are  rational  functions  of  < It  i3  interesting 

op 

to  note  that  if  K(s)  is  a rational  function  of  s,  y y «3  * sn.c>)  is 

r -sT  “ * 

a rational  of  obtained  by  resolving  K into  simple  terms  by  a 

partial  fraction  expansion  and  then  applying  the  summation  used  in  Eq.  6 

to  each  term.  Since  one  goes  throu  h the  same  process  to  relate ^ (s) 

to»«£(s)  as  to  relate  K(s)  to  ^ 2 X(s  + jnXl),  we  shall  refer  to 

n 

1 y Sag.  . 

f K(s  ♦ jn.fl)  as  K (s). 


•00 


71  - -OO 
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2.1  Ei; spouse  of  the  System 

The  above  equivalence  allows  the  immediate  reduction  of  any  cascaded 

combination  of  components  to  the  standard  form  of  a continuous  filter  (K  ) 

a 

followed  xn  succession  by  a modulator  (M),  a discrete  filter  (D),  and  a second 
continuous  filter  (J^) . (In  degenerate  forms  some  of  these  components  may 
ba  missing* ) The  reduction  proceeds  as  follows:  I i*  ^h;i  cascade  contains 

only  one  modulator,  it  is  already  in  standard  form.  If  it  contains  two 
modulators,  between  t*em  there  is  a transition  from  sanpled  signals  to 
continuous  signals  in  preparation  for  a new  sampling.  Sampled  signals 
are  delivered  to  a continuous  filter  to  be  smoothed  and  are  then  remodulated. 
This  smoothing-rewodulating  operation  then  can  be  combined  in  one  equivalent 
discrete  filter  operation.  Cascade  arrangements  with  any  number  of  modulators 
are  treated  by  applying  this  reduction  to  one  filter-modulator  set  at  a time* 

It  is  true  that  ar,v  sampled-daba  system  with  or  without  feedback 
can  be  reduced  to  a standard  form  containing  nothing  but  parallel  branches 
of  the  above  basic  K&bDK^  type.  At  this  point  we  do  not  wish  to  go  into  any 
details  of  these  reductions  to  equivalent  forms.  (See  Section  2.2).  Here 
it  is  sufficient  to  note  that  solution  of  the  w,  system  response  problem 
in  effect  solves  the  general  analysis  problem,  once  suitable  reduced  equivalent 
systems  are  obtained. 

2.11  An  Exact  Calculation  of  Response 
Suppose  an  input  signal  0^(t)  is  applied  to  the  K^KDK^  system  as 
shown  in  Fig.  8a.  We  denote  the  transform  of  the  input  signal  byt€i(3)  and 
apply  similar  notation  for  the  output  signal.  But  the  same  signal  appears 
at  point  A if  we  think  of  ©^(t)  as  being  generated  by  a single  unit 
impulse  applied  at  time  t * 0 to  a filter  whose  transfer  .function 
is  ^(s).  Thus  Fig.  8b  shows  a system  whose  unit  impulse 
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Fig.  8 Seduction  of  a cascade  systs 


response  is  the  required  {^(t).  We  now  recognize  the  combination  of  a 
continuous  signal  filter  followed  by  an  impale?  modulator  to  which  a sampled 
signal  is  applied  — namely,  the  'unit  impulse.  This  combination  is  then 


reduced  as  shown  in  Fig,  8c  to  the  equivalent  discrete  signal  filter, 

~^>  tot . denoted  conveniently  by  the  starred  notation  . Finally, 

we  determine  the  output  transform  by  multiplying  the  input  transform,  1,  by 
the  overall  system  transfer  function . Tha  result  is  21^*  from 

which  the  output  time  function  can  be  determined. 

The  difference  between  this  procedure  and  the  transform  an.ilysis 
of  conventional  systems  lies  in' the  fact  that  here  the  presence  of  the  modulator 
and  the  multiple  frequency  translations  which  it  introduces  require  the 
summation  process  or  in  addition  to  a multiplication  of  component  transfer 

functions.  Another  way  to  look  at  this  situation  is  to  view  the  modulator 
as  a barrier.  The  barrier  can  be  Jumped  by  changing  to  a time-domain  description, 
sampling  the  input  to  the  modulator,  deducing  the  transform  of  tha  samples, 
and  finishing  in  the  frequency  domain.  This  time-domain  viewpoint  often 
auoras  an  easy  way  uo  rina  without  actually  summing. 
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The  features  of  the  basis  response  calculations  are  Illustrated 

in  the  example  of  Fig.  9.  The  system  of  Fig.  5 may  have  been  derived 

by  reducing  a more  general  system  to  the  standard  form.  Its  response  to  a 
step  input  applied  at  a sampling  time  will  now  be  computed. 


Fig.  9 A basic  system  responding  to  a step 


The  signal  applied  to  the  modulator  has  the  transform « bTb  a)* 

After  modulation  the  signal  has  the  transform  £0;  Sa ) givsn  by, 
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Instead  of  summing  directly  we  may  change  momentarily  to  the  time  domain 

to  3ump  the  modulator  barrier.  The  applied  signal  has  the  transform 
1 1 

ld<I  ■ — — - — | — , which  in  partial  fraction  form  indicates  a step 

i a s s ♦ a ’ * r 
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and  an  exponential  component  of  amplitudes  — and  - i • Tlie  bequ«r>.oe  of 

a £ 

samples  from  the  step  (assuming  that  the  sampling  occurs  just  after  a 
discontinuity/  j&  just  a sequence  cr  samples  all  of  amplitude  — with  the 


transform 


t *j  m f 

i a 1 a 


jjL  e“sT  ♦ e"2sT  * .., j 
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The  sequence  of  samples  from  tne  exponential  similarly  has  the  transform 

to,w  \*  m 1 ■)  . -aT  -sT  -2aT  -2sT  . 

\i8fc  1 ♦ e « * e e a...  f oo  • 

x a.  x n 


Rolling  up  these  expressions  into  closed  form  by  inspection  and  adding,  since 
the  modulation  process  is  linear*,  we  find  the  total  modulator  output 

• X* 
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This  answer  is  the  same  as  that  found  previously  by  direct  Summation. 

To  find  the  output  signal  transform,  we  now  treat  the  remainder 
of  the  system  in  conventional  manner  and  multiply  the  remaining  transfer 
functions. 
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x'o  get  an  overall  picture  of  this  response  we  plot  the  pciee 


and  zeros  ofs©s(s)  in  the  s-plane,  as  in  Fig.  10. 
c 


In  general,  several  infinite  sequences  of  poles  and  seres  are 

“*sT 

obtained  because  of  the  factors  of  *€^(s)  which  are  rational  in  e . In 
addition,  the  re  are  factors  rational  in  s which  lead  to  a finite  number  of 
poles  and  zeros.  The  e~  factors  yield  a pattern  of  period  jSl  . 

In  this  diagram  various  components  of  the  output  may  be  associated 
with  the  poles.  The  pole  at  the  origin  signifies  tnat  part  of  the  output 
is  a step  function.  The  two  other  poles  on  the  real,  axis  represent  two 
exponential  transient  components  with  different  time-constants . These  three 
poles  together  form  a central  group  of  relatively  low-frequency  response. 

It  is  for  this  part  of  the  response  that  a sampled-data  system  is  designed. 
The  remaining  polos  represent  that  part  of  the  response  called  the  ripple. 
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Tha  sequence  of  Doles  at  ♦j fi,  *2  j -fd.,  ...  represent  a ripple  component 

persisting  in  time,  the  steady-state  ripple.  The  poles  at  ♦jfl.-a,  ... 

represent  a danped  ripple  component. 

Suppose  that  we  let  ? — y 0 - i.e..  J~L — — -y  oo.  Then  the  sanpled-data 

system  approaches  a continuous  system.  In  Fig.  10  the  s - *jil  lines  recede, 

leaving  the  central  group  of  response  poles  in  complete  command  of  the  s-plane. 

We  look  for  our  basic  response  poles  inside  the  strip  shaded  in  Fig.  10 

between  the  lines  e - *jA  and  s * -j  n. 

A sketch  of  the  pole-zero  diagram  of  the  output  transform,  then, 

immediately  shows  us  the  qualitative  nature  of  the  response.  In  addition 

it  suggests  the  application  of  the  procedure  of  numerical  evaluation  of  the 

pole  residues  to  build  up  the  output  response  from  its  various  natural  modes. 

lowever,  an  annoying  feature  of  such  quantitative  work  is  tie  occurrence  of 

-aT 

Hiixed  expressions  containing  both  s and  e t In  the  Sec.  2.12  we  show 

how  one  can  arrive  at  an  approximate  determination  of  response  which  is  simpler 

than  the  direct  method  above-  The  s -plane  diagram  has  its  major  value 

in  exhibiting  an  overall  qualitative  picture  of  the  response.  A sketch  of 

the  s-plane  diagram  would  be  a preliminary  step  in  quantitative  analysis. 

2.12  An  Approximate  Calculation  of  Response 

The  approximate  method  for  computing  the  response  of  the  K MDIC 

a b 

system  rests  on  the  observation  of  the  previous  section  that  the  response 
consist*  essentially  cf  two  parts,  a "first  order  response"  with  components 
of  low  natural  frequency  and  a ripple  response  of  high-frequency  components. 

Our  approximate  method  will  determiriH  each  of  these  effects  separately 
in  the  following  manner. 
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The  first  order  response  may  be  found  by  computing  the  system 
output  only  at  tne  instants  at  which  the  impulse  modulator  sample  s its 
input*  In  other  words,  we  inquire  what  the  output  is  at  t = 0,  T,  2T,  3T... 
tc.  a transient  input  applied  at  t « 0.  This  procedure  is  approximate  only 
in  the  sense  that  we  ignore  what  the  output  does  between  the  sampling 
instants^  the'values  of  the  output  at  these  instants  are  to  be  determined 
accurately. 

In  order  to  determine  the  output  of  the  X MDK  system,  we  connect 

a d 

to  the  output  of  the  system  a second  impulse  modulator  which  runs  in  phase 
with  the  first  (see  Fig.  ll). 

*>  H--IEI— S3 — ■ 


b)  HHS-0- 


original  system 
modified  system 


c) 


equivalent  of  modified 
s vs  tern 


^o  " HV 


Fig.  11  Reduction  of  second  modulator  from  a cascade. 


Then  the  new  output  is  a sequence  of  impulses  whose  areas  are  identified 
with  the  value  of  the  original  output  at  the  sampling  times.  Since  at 
point  A and  the  new  output  terminal  the  signals  are  of  sampled  nature, 
an  equivalent  discrete  signal  filter  maybe  found.  Its  traisfer  function  is 
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the  tranjform  of  the  sampled  output  when  the  input  transform  is  given,  we 
obtain  immediately*^  * DK^  * whert5  the  sign  in  the  symbol  *9^ 

indicates  that  ve  are  dealing  with  the  .sanoled  output.  This  result  differs 
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^ *x* 

from  the  previous  one  in  that  nnai^  is  rational  function  of  o ° with 

no  factors  rational  in  s like  those  which  occurred  in  the  expression  for . 

o 

In  the  previous  discussion  the  s-plane  wse  us-.d  to  sketch  the  poles  and  zeros 

-sT 

01  the  response.  Here,  however,  it  is  more  useful  to  use  the  e plane* 

-_T 

Furthermore,  since is  rational  in  e , only  a finite  number  of  poles 

-sT 

and  zeros  appear  in  the  e plane.  As  far  as  the  sampled  output  is  concerned, 

then,  only  a finite  number  of  modes  of  vibration  arise.  Their  magnitude 

-sT 

may  be  found  by  either  sample  graphical  measurements  on  the  e plane 

* 

pole-aero  pattern  or  by  breaking  up  analytically  into  partial  fractions, 
each  term  of  which  represents  a single  mode  of  response. 

The  analytical  procedure  is  il3ustrated  in  the  following  example. 

The  diagram  of  the  system  to  be  investigated  appears  in  Fig.  12.  A unit 
step  will  be  applied  to  the  input.  We  renuire  the  output,  at  the  sampling 
instants. 
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Fig*  12  Process  for  finding  sampled  response 
The  complete  output  has  the  transform 
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The  sssplea  output  has  the  transform 


m*  - (|a )*  (i  - e"aT) 


-pT 


- I 

s(s+a)  J 


-sT\ 


a-.-*)  .-,T 


( 2$) 


*fci\ 

itJL-e  y 


(e~sV 


/.  -sT'v>  aT  -sTx 
(1-e  ) (1-e-  e ) 


The  character  of  the  output  is  indicated  by  plotting  the  poles 

and  zeros  ofK*^  in  the  s-piane  (sea  Fig.  13a).  The  transient  frequencies 

of  the  sampled  response  (first  order  system  response ) arc  indicated  by  the 

4*  “*sT 

pole-sero  diagram  of^  in  the  e plane*  (Sv#  Fig.  13b). 
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Fig.. 13  Pole-zero  pattern 
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The  s-plane  diagram  indicates  that  the  overall  response  has  the  following 

components : a linear  runp  component  indicated  by  the  double  pole  at  the  origin} 

possibly  a step  component  (whose  simple  pole  is  obscured  by  the  double  pole)} 

a de saying  exponential  transient  indicated  by  the  pole  on  the  real  axis}. and 

a constant  steady-state  ripple,  indicated  by  the  sequence  of  simple  poles 

on  the  imaginary  axis.  Transient  ripple  is  absent. 

-sT 

The  e plane  diagram  shows  similarly  that  the  sampled  response 

**sT 

indicates  a linear  ran p (double  pole  at  e “ l),  possibly  a step,  and  an 

exponential  transient  (simple  pole  at  e”3*  ■ e ax).  We  can  find  these 

* 

components  by  breaking  up  the  expression  into  a partial  fraction  expansion 
which  exhibits  these  modes.  In  order  to  see  what  to  look  for.  we  have 
recorded  in  Fig,  ili  some  simple  nodes  and  their  transforms. 
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Fig.  lii  Time  functions  and  their  transforms 
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Tha src  have  fouiid  breaks  up  in'.o  Appropriate  fractions  in  the  following 
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Note  that  contrary  to  the  usual  practice  in  suen  expansions  we  look  for  a 

e‘sT  1 

terra  of  the  form  — ^ 2 rather  than  ~-a!~2  * 3ince  the  former  can 

(1-e  ) (1-e”8  ) 

immediately  be  identified  as  the  transform  of  a sampled  ramp  function. 

This  change  does  not  occasion  any  extra  computational  difficulty*  Now 

we  identify  the  first  term  of  the  result  is  a sampled  unit  linear  ramp, 

T 

the  second  terra  as  a step  of  amplitude  - , and  the  third  is  a 

1 T 

decaying  exponential  of  time-constant  — and  Initial  amplitude 


T? * 

a 

For  c6uvenier.ce,  then,  the  sampled  output  may  be  considered  as  arising 

from  the  seap?..lhg  of  the  continuous  signal  0 (t)  * t - - — — m — (l-e~aT), 

8 1— e 

and  this  signal  may  be  taken  as  an  approximation  of  the  actual  response 

9^(t)  where  rippl«*  effect*  have  been  neglected* 

A sketch  of  the  input  to  this  system  9^(t),  the  actual  response 

©Q(t),  the  sampled  response  ©Q*(t),  and  the  envelops  of  th$  sampled  response 

9 (t)  appear  in  Fig.  15. 

6 
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Fig,  15  Plot  of  approximations  to  the  response. 

The  envelope  curve  of  the  sanples  approximates  quite  Kell  the  basic  first 
order  behavior  of  the  system „ We  are  left  with  the  problem  of  estimating 
the  ripple  - of  finding  out  what  happens  between  saTTples.  A glance  at  the 
response  curves  of  Fig-  15  shows  t.hat.  if  y#  were  +o  sample  the  output  every 
T/2  seconds  instead  of  only  every  T seconds,  the  samples  we  wouj.  1 get  would 
show  a component  of  oscillation  of  frequency  -Cl , The  amplitude  would  indicate 
very  wall  the  amount  of  ripple.  As  long  as  the  system  under  consideration 
has  such  a simple  intersample  behavior,  this  ripple  estimation  will  be 
reasonably  accurate r.  In  most  control  systems,  since  the  controlled  output 
member  has  low  pass  filter  characteristics,  the  method  may  be  satisfactorily 
applied. 

Suppose  we  wish  to  double-frequency  sample  the  output  of  some 
system.  A second  impulse  modulator  operating  at  one-half  the  period 
of  the  first  one  is  attached  to  the  output,  and  is  synchronized  with  the 
first  modulator  (see  Fig,  16) . 
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Fig.  16  Block  diagram  ol'  system  for  rippls  calculation.. 

Now  at  point  A of  Fig.  16a  tue  signals  are  sampled^  one  sample  arriving 
every  T seconds.  Nothing  prevents  us,  however,  from  considering  these 
samples  to  be  arriving  on«;  every  T/2  seconds.  We  Just  say  that  every 
other  sample  ties  a value  of  zero.  Then  «+•  point  A and  the  sampled  output, 
the  signals  are  sampled  with  period  T/2  and  an  equivalent  discrete  signal 
filter  exists  given  by 

) (27) 


can  in  most  practical  instants  be  found  directly  flrom 

♦ ^ 

the  used  in  determining  the  first  order  response.  is  a rational 

-sT  _ **  -sT/2 

function  of  e • similarly  is  a rational  function  of  e ' . 

may  be  found  from  merely  by  changing  every  T to  T/2,  provided 

that  the  correct  square  root  is  obtained.  For  instance,  if  somewhere 


-sT 


the 


in  the  expression  for  some  e was  obtained  actually  by  reducing 
more  -general  expression  e 1 where  the  system  constant,  a.  Just  happened 

to  have  the  value  then  Just  replacing  T by  T/2  would  give  the  result 

-aT/2  -(s*a)T/2 

e ' in  the  case  of  the  reduced  expression  and  e in  the  second. 

Put  the  former  result  is  not  now  correct,  3ince  e”a^^  » e”*^,  which  is  not 

equal  to  1 . Thus  the  simple  replacement  of  T by  T/2  may  be  performed  only 


'.1 L.  r»  rtrtA 

n&yjrh  wce.c 


ir  such  reductions  have  not  been  made  either  because  the  original  expression 

was  retained  or  because  the  original  expression  could  not  be  reduced  at  *11, 

Of  course,  we  way  always  return  to  the  full  expression  so  as  to 

be  sure  what  to  do,  but  bi*e  eight  be  saved  if  «■«*  could  tell  at  a glance, 

given  X^  , whether  X^  could  be  directly  obtained  from  it.  This  simple 

r-t-  \T/o 

substitution  procedure  is  possible  if  the  e *i  expressions  resulting 
froa  the  presence  of  nodes  of  oscillation  at  the  frequencies  s » -e^  cannot 
be  simplified.  The  expressions  cannot  be  simplified  if  these  frequencies 
(indicated  by  polee  of  in  the  s-plane)  lie  in  the  strip  froa  s » -j-fl. 
to  b *»  * j Xl . If  any  poles  lie  on  the  edges  ox  this  strip  or  outside  it, 
the  corresponding  expression  for  5^  contains  factors  of  the  type  e ”J> 
which  if  removed  throw  away  essential  phase  information  needed  if  7/2  is 
substituted  for  T. 

As  an  example,  consider  two  cases.  In  the  first  X^1^  - — ^ • 


in  the  eecond  t 


«)  - 


(a  is  real).  The  pole-xero  diagrams  of  these 


X^  are  shown  in  Fig.  17d.  We  compute  x£^*,  X^‘c^**. 


(1) 


Fig.  17  Pole-xero  diagrams  for  X^ 
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Note  the  difference  in  results.  K.  ' is  an  entirely  different 

o 

expression  that  the  final  result  for  while  may  be  obtained 

(2)* 

from  by  Just  changing  T to  T/2  without  going  back  to  an  earlier 

stage  of  the  computation.  There  is  in  fact  in  the  latter  instance  no  earlier 

stage,  since  no  reduction  of  the  original  summation  is  possible.  Note  how 

this  behavior  corresponds  with  the  pole  positions  of  the  in  the  s-plane. 

Once  has  been  determined  and  the  T/2  period  sampled  output 

transform  0^^  - has  been  found,  the  poles  and  zeros  of 

9q  may  be  located  in  the  e plane.  (Note  that  use  ha?  been  made  so 
-sT  -sT/2 

far  of  the  s,  e -plane  diagrams.  Th*  c -p1  cr.3  is  a new  one.)  If 

V can  be  obtained  from  by  simple  replacement  of  T by  T/2,  the 

-*bT/2 

poles  and  zeros  of  in  the  e ' -plane  are  the  same  as  the  poles 

* -sT 

and  zeros  of  in  the  e -ulane.  These  latter  have  already-  been  found 

during  the  determination  of  the  first-or^er  response.  The  polp?  and  zeros 
.*  -sT/2 

of  (9jK  ) D in  the  e -plane  may  always  be  found  by  taking  both  square 
roots  of  each  pole  and  zero  of  (©^K^)  D in  the  a -plane.  For  example, 
if  (0,5^) ”15  has  a zero  at  9%.  60°  in  the  e 3 1 -plane,  then  In  the 

-cT / ^ O "O 

e “ -plane  there  are  corresponding  zeros  at  3 2$.  30  and  3 2^.210  . 


*!* 


Therefore  the  pole  and  zero  pattern  provided  by  the  (©.,&_)  D part  of  ©q 
"ST/2 

in  the  e -plane  is  symmetric  about  both  real  and  imaginary  axes,  since 

-sT 

the  pattern  in  the  e -plane  is  already  symmetric  about  the  real  axis. 

There  will  arise  in  the  above  manner,  a finite  number  of  poles 
-sT/2 

in  the  e ' -plane,  each  one  corresponding  to  a cojjponeni  of  the  double -sampled 


V M h4  W 


r?e  can,  if  we  like,  regard  this  pattern  as  indicating  a higher 
approximation  to  the  actual  (unsampled)  output  respond®.  However,  the 


si 


first  order  behavior  has  already  been  obtained  through  the  e -plane 


diagram,  arid  if  the  ripple  is  small,  then  this  higher  approximation  will 
not  lead  to  a significant  change  in  the  magnitude*  of  the  first-order  response 
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purpose  of  estimating  the  ripple.  This  ripple  magnitude  is  associated 

.ct/ 2 

with  the  i^sidues  of  poles  in  the  e ' plane  not  already  identified  with 
first-order  response)  components.  Such  residues  may  be  found  analytically 
or  by  graphical  measurement. 

The  above  discussion  is  now  illustrated  by  finding  the  ripple 

« 

in.  the  system  of  Fig.  9»  whose  first  order  response  has  been  found  and 
sketched  in  Fig.  12,  First  note  the  ^ -j  has  poles  within  the 

strip  s - to  s - Therefore  is  found  directly  from  K,  *, 

. , , . _ _ . ,,  -sT/2  plant  is  the  same  as  that  of 

and  the  pattern  of  F in  the  e 1 

* -sT  b ** 

K.  in  the  e plane . ^ may  be  found  immediately. 
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A sketch  of  the  *©^  pattern  in  the  o plane  is  shown  in  Fig.  18. 


(■©k  pattern 

Fig.  18.  Pole-zero  pattern  of  f©^**  in  e”S1/ ^-plaiie. 
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There  ore  four  possible  output  components  -indicated  by  the  poles  in  Fig.  1 

& -ST 

Compare  this  pattern  with  that  of  ^ in  the  e plane  in  Fip - 13b,  The 

o 

at  e = 1 and  e * 0-ai/2  in  Fig.  18  correspond  to  those  in 


pOiOS 

Fig,.  13b  for  which  5 =1  and  e.  - e = and  indicate  so  before  a ramp, 

a otep,  and  a decaying  exponential  component  in  the  output = However,  there 

-»T/p  S-* 

is  now  an  extra  pole  at  e “ ' « -1  in  the  ^ pattern.  This  corresponds 


to  a term 


1+e 


A 

in  the  partial  fraction  expansion  of  and  represents 


^f/2 


(refer  to  Fig.  ll)  a sampled  cosine  oscillation  of  amplitude  A and  frequency 

_TL  . We  note  from  the  s -plane  diagram  of  ^ (Fig.  10a)  that  there  are  no 

transient  ripple  components  and  are  not  surprised  that  we  obtained  only 

the  one  ripple  pole  . Evaluating  A we  get 

-sT 

A < 


(l*e  ‘ ) ,»0* 


-»T 

-T(l-e  g ) 
-aT 

U(l*e  T) 


**. 


-sT 

. F 


-1 


so  that  the  ripple  may  be  considered  as  arising  from  the  samples  of  a wave 
T (] r~\ 

- i cos  Jit,  and  this  wave  may  be  added  to  the  first-order  response 

4 (i+e'aT/2^ 

approximation  to  obtain  a continuous  wave  which  fits  very  well  to  the  actual 
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The  step  component  of  the  first-order  response  should  be  modified,  however, 
either  by  redetermining  it  from  the  Ksd  transform,  or  by  merely  matching 

U 

the  ripple  to  what  is  already  known  to  be  the  value  of  the  samples  at 
t * 0,  T,  2T  ...  • Thus  to  the  first  order  wave  one  should  really  add 


l (i-e-»T/:j  n+\ 


(l+e 


•a?  r? 

— 


in  ni  t - 0,  T,  21  ...  the  samples 


of  the  new  curve  have  the  same  value  as  before.  The  second-order  approximation 
includes  the  ripple  effsct.  The  envelope  curve  of  the  double  frequency 


A ^ 4-Wa  4 0 < 


«.<*>  “ 1 • r^sr. r 


x-e 


( -?] 
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U / -a? ' 
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From  this  erpreasion  we  can  find  how  the  ratio  of  the  step  lag  in  response 
to  the  ripple  amplitude  behaves  with  a change  in  i , the  system  time-constant, 
as  well  as  noting  their  absolute  values.  This  method  of  approximation 
to  the  response  of  the  sampled-data  system,  then,  readily  provides  enough 
information  for  design.  It  offers  means  for  dete mining  just  how  system 
parameters  should  be  set  for  optimum  response  characteristics. 

Summary 

The  preceding  section  discussed  the  analysis  of  the  response  of 
a basic  cascade  sampled-data  system.  It  was  mentioned  that  any  sampled-data 
system  of  arbitrary  complexity  could  be  reduced  to  an  equivalent  system 
containing  at  most  a number  of  parallel  branches  of  the  standard  cascaded 
sequence  of  elements,  continuous  signax  i liter,  impulse  modulator,  discrete 
signal  filter,  and  second  continuous  signal  filter.  Therefore,  the  analysis 
of  the  response  of  -the  basic  sequence  is  sufficient  to  allow  analysis 

of  any  system.  The  response  of  the  system  was  then  calculated  exactly. 

Finally,  methods  of  approximation  to  this  response  were  given  which  had  a 
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practicsl  value  in  being  easily  applied  to  obtain  an  overall  picture  of  the 
main  features  of  the  response  and  their  dependence  on  system  parameters. 

It  remains  to  demonstrate  in  Sec.  2.2  how  the  general  sampled-data  system 
is  reduced  to  a parallel  set  of  standard  branches. 


2.2  Reduction  of  the  General  Sampled  Data  System 

This  section  is  devoted  to  a method  by  which  an  equivalent  network 
is  obtained  for  a general  system  whose  response  is  easy  to  compute  by 


v.rn  4 


2.21  Symbolism 

To  facilitate  discussion  of  the  points  to  follow  it  is  necessary 
tr  introduce  a shorthand  for  dealing  with  sampled-data  systems.  The  reader  is 
probably  familiar  with  the  way  in  which  the  use  of  transfer  functions  in 


the  t ran « f o mat-ion  of  conventional  linear  networks  into  various  circuit- 


equivalents  leads  to  a compactness  of  notation  and  an  appreciation  of  the 
circuit  behavior  which  cannot  be  equalled  by  the  use  of  the  circuit  differential 
equations  themselves. 

Unfortunately  sampled-data  systems  cannot  be  entirely  described 
by  3uch  transfer  functions.  It  has  been  noted  in  previous  sections  that 
the  discrete  and  continuous  filters  may  be  so  described  but  that  the  sampling 
process,  being  essentially  a modulation,  cannot  be  represented  by  « simple 
transfer  function  (at  least  not  by  one  of  the  ordinary  time-invariant  type). 

One  easy  way  to  appreciate  this  fact  is  to  note  that  the  wave  shape  of  the 
output  of  a sampler  (impulse  modulator)  depends  not  only  on  the  wave  shape 
of  the  input,  as  it  would  if  it  were  des crib able  by  an  ordinary  transfer 
function,  but  in  addition  depends  on  the  relative  timing  of  the  input  and 


sampling. 
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However,  the  sampling  process  is  linear,  the  output  for  a sum  of 
inputs  being  equal  to  the  sura  of  the  responses  from  the  inputs  applied  one  at 
a time.  This  superposition  propei-ty  is  an  essential  feature  of  the  methods 
to  follow.  Because  of  this  linearity  of  the  separate  components  of  « sampled-data 
system  (l.e.,  impulse  modulator,  digital  filter,  continuous  filter),  the  system 
itself  is  linear  in  its  entirety.  The  store  of  superposition  methods, 
suitably  extended,  is  therefore  available  for  the  analysis  of  sampled-data 
systems » 

Since  a sampled-data  system  is  linear  but  not  entirely  uescribable 
by  transfer  functions,  the  more  general  concept  of  linear  ' erations  may  be 
introduced.  Every  component  of  a sampled-data  system  operates  linearly  on 
its  input  to  produce  its  output.  Such  an  operation  may  be  indicated  by  a 
juxtaposition  of  input  symbol  and  operational  symbol  as  indicated  in  Fig.  20, 


Linear 

component 


e. 

3 

1 

V J 

L 

_ v° 

1 

✓ 

Fig.  20.  Representation  of  a general  linear  component. 

Note  that  the  order  in  which  the  symbols  appear  is  significant.  0^L  means 
that  L operates  on  This  order  is  chosen  arbitrarily  to  provide  a certain 

isomorphism  of  the  symbolism  with  the  way  in  which  electrical  engineers 
usually  draw  diagrams  — i.e.,  signal  flow  from  left  to  right.  Thus,  consider 
the  cascade  of  linear  components  in  Fig.  21  and  the  corresponding  operational 


expression 
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Fig.  51s  Severs!  linear  elements  in  cascade. 

The  overall  operation  of  a cascade  of  separate  operations  is  indicated 
by  the  symbols  in  the  order  in  which  the  separate  operations  are 

performed.  The  fact  that  the  operations  are  linear  does  not  imply  that  they 
are  commutable,  that  is,  that  •*  It  in  known,  however,  that  if 

L,  and  L„  refer  to  operations  performed  by  continuous  filters  they  are 

ju  C 

commutable,  the  order  of  cascading  ordinary  constant  coefficient  linear 
networks  being  immaterial.  The  impulse  modulating  operation,  however, 
obviously  does  not  commute  with  a continuous  filtering  operation,  since 
the  output  must  be  in  the  form  of  samples  in  the  first  instance  but  in  the 
form  of  a continuous  wave  in  the  second.  It  must  also  be  noted  that  the 
filter  operations  also  cannot  be  commuted  if  there  is  a sampling  operation 
between  them. 

Because  of  the  linearity  of  the  systems  under  consideration,  the 
operation  performed  by  a paralleled  set  of  components  may  be  indicated  by 
the  sum  of  the  separate  operations.  Refer  to  Fig.  22. 
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Fig,  22,  Linear  elements  in  parallel. 

We  now  allocate  a specific  symbol  to  the  operations  performed 
in  a sample d-data  systems 

M,  operation  of  sampling  or  impulse  modulation; 

K,  operation  of  continuous  filtering; 

D,  operation  of  discrete  filtering; 

T,  overall  operation  (transmission). 

Because  cf  the  sampled  signals  which  are  involved,  the  digital 

i 

filtering  D is  always  preceded  by  either  another  digital  filtering  D or 
impulse  modulation  M. 

An  important  restriction  must  be  made  on  the  systems  to  be  considered. 
It  will  be  assumed  that  all  impulse  modulators  are  running  at  the  same 
frequency  and  in  the  same  phase  (or  at  least  very  nearly  in  the  same  phase). 
Furthermore,  it  will  be  assumed  that  no  questions  arise  as  to  whether 
various  impulse  modulators  take  their  samples  just  before  or  just  after 
input-wave  discontinuities  produced  by  previous  sampling  operations. 

The  restriction  on  frequency  and  phase  is  essential;  the  restriction 
on  the  samoling  of  discontinuities  is  not  essential  and  is  introduced 

s /i_  __  _ __  x — j j.  u . *i  _ _ t j j _ * tn  ^ — . 

vnjLjr  tuiiiauciitc  j.i:  i.jjc  iuhum.ui^  uj.scussj.on,  me  neg-Lecuea  euecx.3 

are  easi?.y  taken  into  account;  how  this  can  be  done  is  best  described  later 
after  the  basic  ideas  have  been  presented.  Ei cause  of  the  above  stimulations 
we  can  make  indiscrirainant  use  of  the  symbol  M for  every  sampling  operation 
in  the  sampleu-data  system. 
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Vary  often  a cascade  operation  occurs  in  which  more  than  one 
sampling  appears;  see  Fig.  23 » 


sample  signal? 


T - MDKM 
- M(DI)" 
* MDF* 


Fig.  23*  Reduction  removing  all  but  the  first  modulator. 

In  this  ease  as  discussed  in  See.  ?.  the  system  ary  fcs  simplified  to  on 
equivalent  system  of  a modulator  followed  by  some  discrete  filter  given 


by 


ij  ^ D(  a+ jn  ■/!)£(  a*  jnSl ) 

n— 00 


~n*“oo 


(36) 


Deq(s)  * D(s)K'*(s),  where  ths  * refers  to  the  indicated 

summation.  Use  has  been  made  of  the 
fact  that  D(s)  is  periodic. 

We  use  the  operational  symbols  in  the  same  way  as  the  transfer  function 
symbols  above.  That  is 


Deq  - (DK>*  - DK* 

so  that 

MDKM  - MDK* 


(37) 

(38) 
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2.22  Reduction  of  Ssrles-Farallfll  Sampled-Data  Network. 

The  operational  symbolism  introduced  above  r»ay  be  used  in  the 
finding  oi  equivalent  circuits.  This  simplest  application  of  these  methods 
is  the  reduction  of  a fieries-parallsl  system  to  simpler  form. 

From  the  linearity  of  a saapled-data  system  it  follows  that  the 
distributive  law  of  "multiplication"  holds  with  respect  to  the  operations 
if  correct  order  of  the  factors  is  maintained.  Refer  t.v  Fig.  2lj. 


et  **ili 

- « 


Reduction  of  a series-parallel  network. 
k this  point  it  is  well  to  take  stock  of  the  relationships  involved 
in  operator  manipulations.  Table  I gives  a .summary. . 
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Summary  of  C'perator  Manipulations 


Conventions  in  Notation 


L •»  anv  sampled-data  system  operation 
M - impulse  modulation,  sampling 
E - continuous  signal  filtering 

D ■ discrete  signal  filtering 

# 

X ■ the  derived  discrete  filter  equivalent  operation  for  the 

operation  K followed  by  K when  the  input  signals  are  discrete* 

T ■ the  symbol  for  some  overall  combination  of  operations  which 
relates  an  input  signal  to  a?',  output. 

A « a simple  gain  or  attenuation 

L.L-L,  - operations  L. , 7j  t L,  performed  successively  in  the  order 
123  1,2,3.  * 2 3 

L, ♦L0+L,1  » operations  performed  separately.  results  added. 

Symbolic  Relations 

a)  Commutability 

ME  + KM  - K2Ei 

DK  / KD  D,D.  - DJX 

x 2T  ± 

MD  4 EM,  unless  the  input  signals  are  discrete. 

All  symbols  commute  with  the  gain  operator.  A.. 

b)  Contoination,  factoring 

WV  “ L1L2*  L1L3 

VL3)L1  ‘ W Vl 


^Vl2)<W  - WWWVu 


c)  Cascade  Equivalents 


(DE)”  - HtT 


>_  * 


'Vz'  * h 


MDEM  - MDK 


Through  use  ef  these  symbols,  any  serles-parallel  samoled-data 
systen  operation  may  be  written  in  symbolic  form.  Reduction  to  a simpler 


equivalent  operation,  if  such  exists,  may  be  effected  by  using  the  relations 
mentioned  abovo.  Fig.  25  gives  an  example  of  such  a reduction.  For  simplicity 
the  blocks  of  the  block  diagram  are  omitted.  Summing  points  are  indicated 
by  nodes,  the  output  on  each  of  the  lines  from  a node  being  the  sum  of  all 
the  inputs  to  that  nCvC  * ^pc  rsoicixS  are  written  above  corresponding 


branches  between  node*. 


(a)  Original  system 


t • ( w. A YiT  fir  .hfr,  r uj.t\  r Nr  1 



- fM-frAl  Pit  If  xK"  MT1  X KIT  <*V  fflf  K.  X . 

...  L-r-2  -1— i-j-5  -1—  2TU"5J 

- (m*a)  [k1k2*k1md1k3*s5+k1md2iuk5] 

or  - HF1K2*MK1*D1K^F.^MK1*D2EliK5 


♦ AK1K2+AK1MD1K3  ^♦AK1KD2^Kcr 


0 


(b)  Equivalent  system  in  standard  form. 

Fig.-  25.  Reduction  of  a complicated  system  to  standard  form. 
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The  results  of  this  reduction  wav  not  at  first  glance  appeal' 
simpler  than  the  original  expression.  In  terms  of  economy  of  symbols 
this  observation  is  correct.  However,  the  object  of  practical  interest 
is  to  t-e  able  to  compute  the  output  from  such  a system  when  the  input  is 
given.  The  new  form  permits  this  computation)  the  old  form  doas  not,  unless 
one  simultaneously  during  the  computation  goes  through  the  various  steps 
mentally  by  which  the  equivalent  was  obtained. 

Notice  from  the  example  that  a series-parallel  network  can  always 
be  reduced  to  a parallel  3et  of  operations  of  the  form  K KISS^ . (The  operation 
E is  a special  case  of  this  form,  when  M and  D are  absent.)  To  compute 
the  overall  output,  it  is  only  necessary  to  sum  each  of  the  outputs  from 
each  operation 

The  computation  of  the  response  to  the  operation  has  already 

been  described  in  Sec.  2.11,  where  it  was  fihovn  Kb  at.  given  the  transfer 
functions  V s),  K^(a),  D(s)  and  the  signal  transformers),  the  output  is. 


*«(a)  - ^(s+jnir).)a1(s'*'jn/l) 

”n»- 


) D(8)K?rs) 


D(s)K2(s). 
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Every  term  in  the  reduced  series-parallel  equivalent  can  be  handled 
in  this  way.  Note  that  if  two  or  more  terms  of  the  reduced  equivalent  have 
the  sane  operations  before  (or  after)  the  sampling,  a rsctoring  is  possible 
which  still  preserves  the  simplicity.  For  instance  suppose 

T - 

“ ^(D.VV;*),  (U0) 

so  that  the  output  raay  be  computed  as  follows, 

r„\  _ rL./_w  /_n  * /.%»  (a\  i n fJT!  i'Im 

Hr  \a/  5 I \o/nA^/  T V4J*/ 

° L1  1 J L!  ! ‘ 3 J 

which  is  seen  to  be  the  same  result  that  would  have  been  obtained  from  finding 
the  responses  separately. 

2.3  Reduction  of  Feedback  Sampled-Data  Networks 


The  reduction  of  series-parallel  networks  is  very  straightforward 
since  the  overall  operation,  T,  say  immediately  be  written  by  inspection  and 
subsequently  reduced.  It  is  not  so  obvious,  however,  how  to  write  down  an 
expression  for  the  overall  operation  performed  by  a feedback  system,  since 
there  are  constraint  relationships  among  the  various  signals in  the  system. 
The  purpose  of  this  section  is  to  show  how  any  feedback  system  can  be  re- 
duced t.o  a series -parallel  system  whose  output  can  be  directly  evaluated 
in  the  manner  previously  described. 

The  entire  procedure  hinges  on  two  conditions  - the  linearity 
cf  the  system  in  general,  and  the  possibility  in  particular  of  finding 
a directly-connected  circuit  which  performs  the  same  operation  as  a feedback 
loop  containing  a series-parallel  combination  of  sampled-data  system 


components 
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The  remainder  of  the  discussion  makes  use  of  flow-graph  methods. 
The  translation  from  block-diagram  to  flow-graph  language  in  shown  in 
Fig.  26. 

Block  Diagrams  Flow  Graphs 
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Fig.  26  Block  diagrams  and  their  corresponding  flow  graphs. 


The  flow  graph  is  simply  a convenient  way  to  draw  block  diagrams. 


Actually  every  operation  which  involves  manipulation  of  p flow  graph 


could  be  applied  equally  well  to  a block  d-ingram.  The  essential  idea  here 


is  that  block-diagram  or  flov-graph  techniques  can  be  applied  equally 
well  to  servo  systems,  feedback  amplifiers,  linear  computing  systems, 


TV  ..  . TA  AAA 

xtfsuOFu  n.—ccrt 


-M- 

yi  cry  ether  system  comprised  of  lines  •*  operations . une  can  snow 

interdependences  of  signals  with  flow  graphs  for  non-linear  systems, 

hut  the  reduction  manipulations  do  not  apply  because  the  reductions 

7 

are  based  cn  superposition. 

The  problem  in  reducing  a flow  graph  is  to  reduce  the  number 
of  nodes  to  a minimal  number  and  to  reduce  the  number  of  branches  to  a 
minimal  number.  For  a conventional  linesr  system  a fully  reduced  flow 
graph  is  one  branch  connecting  the  input  node  and  the  output  node* 
Reductions  are  of  three  kinds:  (l)  A node  between  two  cascaded  elements 

can  be  removed^  (2)  parallel  branches  between  the  same  pair  of  nodes 
can  be  replaced  by  a single  branch^  (3)  reductions  cf  self-loops  and 
their  nedss  can  be?  made.  when  one  has  made  as  many  reductions  of 
types  1 and  2 as  possible,  there  will  be  a number  of  residual  nodes 
and  two  kinds  of  loops.  Self-loops  will  be  made  of  une  branch  and 
will  start  and  end  on  the  .same  node,  a?;  in  Fig.  28a  = System  loops  will 
involve  several  branches  and  several  nodes.  By  removing  self-loops 
and  their  nodes,  one  can  make  successive  reductions  until  the  flow 
graph  is  completely  reduced  and  has  no  loops  of  either  type.  By  this 
procedure,  all  loops  are  reduced  to  self-loops  and  then  removed.  In 
any  specific  problem,  the  reduction  can  proceed  in  several  fashions, 
but  it  is  well  to  use  the  one  which  involves  least  work. 


Nason,  S.J.,  Some  Properties  of  Signal  Flow  Graphs  (to  be  published 
in  IRE  Proceedings). 
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CLeps  i and  '?  can  usually  be  done  at  once  by  picking  the  residual 
nodes  and  then  evaluating  the  transmissions  between  them  by  inspection. 

For  example,  consider  the  feedback  network  of  Fig.  27. 


Fig.  27.  A block  diagram  cuiu  -Lt/3  flow  graph  which  will  be  reduced. 

Because  of  linearity  we  can  pick  any  nodes,  say  the  circled  ones  A,B,C,  and, 
by  writing  the  relationships  imposed  on  the  signals  which  leave  one  node 
and  arriv-  at.  snothe  •,  the  remaining  nodes  being  left  out  of  the  picture 
for  that  determination,  can  draw  an  equivalent  circuit  containing  only 


the  circled  nodes. 

In 

the  example: 

A to  A 

, 0 

B to  A 

, 0 

C 

to 

A 

, 0 

A to  B 

, 

B to  B 

, (k2k3*i)kukx 

c 

to 

B 

, 0 

A to  C 

, 0 

B to  C 

tr 

c 

to 

C 

, 0 

so  that  the  reduction  to  Fig.  28a  may  be  made  . 


Fig.  28.  The  reduced  graph  for  Fig.  27? 
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Since  the  operations  indicated  are  continuous  signal  filtering.,  we  may  look 
on  the  operation  K as  the  transfer  function  K and  apply  a well-known  result 
that  the  loop  of  Fig.  23  is  equivalent  to  a directly  connected  element  whC3e 
transfer  function  is  i (r  g^el)K  overa13  response  of  the  network 

is  then  given  by  T ■ Vg , See  Fig.  28b. 

Two  essential  features  of  this  process  should  be  noticed.  One 
is  that  the  nodes  to  be  retained  - residual  nodes  - are  chosen  sc  that  the 
operations  on  signals  flowing  among  these  nodes  can  be  written  directly 
in  series -parallel  form.  In  other  words,  the  nodes  are  so  chosen  that  their 
deletion  from  the  diagram  interrupts  all  the  feedback  paths.  Of  course, 
the  input,  and  output  nodes  are  also  retained.  The  other  is  that  an  equivalent 
direct  operation  must  be  found  for  the  self-loop  feedback  loops  which  arise 
in  the  residual  node  graph.  In  conventional,  continuous-signal  systems, 
this  can  always  be  done  by  manipulating  transfer  functions. 

In  sampled-data  systems,  the  picture  is  slightly  different. 

Ir.  the  example  up  to  the  point  where  the  feedback  loop  was  removed  from 
the  residual  graph  there  need  have  been  no  distinction  between  the  treatment 
of  conventional  and  sampled-data  systems.  The  various  K's  could  just  as 
easily  have  been  any  sampled-data  system  operators,  provided  proper  attention 
was  paid  to  the  ordering  of  the  symbols.  However,  since  there  is  no  conventional 
transfer  function  which  describes  the  general  sampled-data  system  operator, 
the  method  applied  to  remove  the  loop  in  the  conventional  system  Cannot 
he  used  in  sampled-data  systems  unless  the  loop  operations  are  describable 
by  transfer  functions } that  is,  if  the  loop  transmission  may  be  described 
as  a discrete  signal  filter  or  as  a continuous-signal  filter.  In  this 
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latter  Instance  the  reduction  can  proceed  in  the  usual  manner  with  the 


interpretation  that  the  transfer  function,  |,~‘p^p.T»  has  a3  its  counterpart 
an  operator  T^-  « — . Note,  however,  that  when  the  loop  contains  a 

sampling  operation  without  a sampled  end  signal  (e.g.,  T - KM)  there  is 
no  transfer  function  • and  the  operation  makes  no  sense  in 

the  above  respect. 

V.'e  are  left,  therefore,  with  the  task  of  finding  some  direct 
operation  with  which  to  replace  the  loop  operations  occurring  in  the 
residual  graph  of  a sampled-clata  system  in  the  general  case.  We  begin 
by  considering  the  simple  case  shown  in  Fig.  2?  for  the  loop  operation 
KjHDKg. 


Ill 


Fig.  ?9.  Handling  a self -loop  with  an  impulse  modulator. 

Instead  of  accepting  the  form  of  the  graph  as  shown  in  Fig.  29a 
we  make  the  following  modification.  As  residual  nodes  of  a new  graph  we 
pick  the  input,  the  output,  and  the  nods  Just  following  the  sampling 
operation  in  the  original  graph.  This  set  is  evidently  satisfactory,  since 
thm  feedback  has  been  removed  through  this  la3t  choice.  The  new  residual 
graph  is  shown  in  Fig.  29c. 
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Now  the  signal  o entering  the  residual  loop  node  are  sampled. 

Also  the  last  operation  which  appears  in  the  loop  is  a sampling.  Therefore* 
an  equivalent  discrete  filter  may  be  found  by  a process  previously  described. 

for  sampled  inputs.  (U2) 

Now  the  loop  operations  may  be  described  by  transfer  functions,  and  sense 

may  be  given  to  the  operation  The  network  has  thus  been  reduced 

l-DdjKg) 

a series-parallel  form  which  is  without  further  transformation  in  the 

standard  form  of  a parallel  set  of  K MT)k  branches.  The  response  is  directly 

a b 

calculable  by  the  methods  of  the  preceding  section,  once  the  input  is  prescribed. 

An  elaboration  of  the  above  procedure  permits  the  reduction  of  any 
feedback  network  to  .series-parallel  form.  The  steps  to  be  followed  ares 

1.  Choose  residual  nodes  in  the  following  way.  Use  the  input  and 
output  nodes  ss  residual  nodes.  Break  up  as  many  feedback  paths  as  possible 
by  picking  residual  nodes  at.  points  in  the  graph  at  which  the  signal  is 

of  a purely  sampled  nature.  (The  insertion  of  branches  of  unity  transmission 
may  be  necessary.)  If  this  procedure  does  not  result  in  the  interrupting 
of  all  feedback,  it  must  be  because  there  are  some  feedback  paths  left 
along  which  all  the  components  are  continuous  filters.  Pick  enough  extra 
residual  nodes  to  interrupt  this  feedback.  Note  that  there  are  three 
classes  of  residual  nodess  a)  input-output*  b)  nodes  at  which  the  signals 
are  purely  sampled,  and  c)  nodes  whose  self-loop  transmissions  can  be 
described  by  a product  of  continuous-signal  filter  transfer  functions, 

2.  Simplify  the  graph  by  inserting  the  equivalent  straight 
through  transmission  for  the  loop  transmission  of  one  of  the  nodes  of 
class  c.  Re-evaluate  the  transmissions  ameng  all  remaining  nodes.  The 
resulting  graph  has  one  less  node  of  class  c but  has  the  original  number  of 
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class  b nodes.  Nov  repeat  this  reduction  until  the  class  c nodes  are  removed. 
Only  class  b ana  a.  nodes  remain. 

3.  Simplify  the  graph  further  by  inserting  the  equivalent 
straight-through  transmission  for  the  loop  transmission  of  one  of  the  class  b 
nodes.  This  step  is  possible  since  the  loop  transmission  is  of  the  form 

c?srating  on  purely  sampled  signals.  Therefore, 
an  equivalent  loop-transmission  transfer  function  eAists  given  by 
(K1  Wr-3  ...),  and  the  equivalent  straight-through  transfer  function 

is  5^ . Now  re-evaluate  the  transmissions  among  all  remaining 

l-o^  d3k2+  ...) 

nodee.  The  resulting  graph  has  one  less  node  of  class  b.  Repeat  this  process 
and  remove  all  class  b nodes.  The  only  residual  nodes  left  are  the  input 
and  output;  the  feedback  has  been  removed  and  the  graph  io  now  in  series-parallel 
form.  The  canonic  KMDF  form  can  now  be  found  and  the  response  to  an  input 
computed. 

Note  that  in  the  above  process  it  is  essential  to  start  removing 
class  c nodes  first  ar.d  then  proceed  to  class  b residual  nodes.  Otherwise, 
if  a class  b node  is  removed  first,  the  loop  transjuiasions  in  the  reduced 
graph  may  be  such  that  some  of  original  class  c nodes  lose  their  charaoter. 
However,  if  a class  c node  is  removed  the  signals  at  class  b nodes  remain 
sampled  and  class  b nodes  remain  as  such  in  the  reduced  graph. 
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As  an  Illustration  of  the  Above  method's  consider  the  saispled-data 
feedback  system  of  Fig.  30«: 
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Fig.  30.  Reduction  of  & flow  graph  for  a complicated  a ampler) -data  system. 
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The  system  in  flow  gr-ijin  form  5 s shown  in  Fig.  jub.  An  extra  branch  of  unity 
transmission  is  inserted  as  in  30c  so  that  a node  b may  be  chosen  which 
interrupts  feedback  and  to  which  the  signals  are  purely  sampled.  There  remains 
a feedback  path  which  cannot  be  broken  by  such  a clas3  b node.  A class  c 
residual  node  is  assigned  which  completes  the  break-up  of  feedback.  In 
Fig.  30c  these  residual  nodes  are  circled  and  their  class  is  noted.  Fig.  30d 
is  the  reduced  graph  showing  only  these  nodes.  The  equivalent,  straight-through 
connection  is  made  for  node  c in  Fig.  30ef  and  the  graph  is  reduced  to 
that  of  30f.  The  only  class  c node  has  been  removed^  now  the  class  b 
node  is  removed  in  the  same  way  and  the  graph  simplified  in  steps  g,h-. 

The  final  result  is  a network  in  standard  form  which  permits  response 
computations  according  to  the  methods  of  Section  2.22. 

• 2.31  A Samplsd-Uata  •Ser YOmecndiiii.Bni 

One  sampled-data  feedback  system  of  particular  interest  is  the 
simple  servo  system  shown  In  Fig.  31= 


Load 

Disturbance  0 


Fig.  31 • Block  diagram  of  a single-loop  servo  system  of  particular 
interest. 

In  this  diagram  the  impulse  modulator  is  followed  by  a filter  with  transfer 
function  K^(s)  which  smooths  out  the  discrete  error  signals  from  the  modulator 
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and  delivers  this  error  to  the  controlled  element  described  by  the  transfer 
function  K^Cs).  Input,  output,  and  load-di3turbance  time  functions  and 
their  Laplace  transforms  are  noted  on  the  diagram. 

Two  response  characteristics  of  this  system  are  of  interest  - the 
output,  resulting  from  a given  input  and  that  resulting  from  a given  load 
disturbance.  In  order  to  make  possible  the  calculation  of  these  responses 
we  transform  the  given  circuit  into  one  without  feedback.  Fig.  3?  shows 
the  application  of  the  general  procedure  for  reduction  given  above  cxf!  minating 
in  networks  of  standard  form. 


1 

1 M i+(K.;K2)*  K.Kp  0 


(a)  Equivalent  System  for 
Input-Output  Response 


Load  Disturbance 


(b;  Equivalent  System  for  Load 
Disturbance-Output  Response 


■fig.  32.  Reduction  of  the  flow  graph  of  the  single-loop  servo  system. 
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As  an  alternative  to  this  procedure  we  may  simply  juggle  the  binds 
and  signals  in  the  system  block  diagram,  making  full  use  of  linearity# 

Although  this  procedure  leads  in  more  general  networks  to  very  cumbersome 
manipulations,  the  above  system  is  simple  enough  to  indicate  its  utility. 

It  should  be  understood  that  nothing  new  is  being  added.  The  following 
method  is  less  sophisticated  but  mechanically  more  involved  than  the  procedure 

for  reduction  already  given  to  which  it  is  basically  equivalent, 
a)  a)  .“'.head  msturoance 


fig.  33..  Block  diagram  manipulations  Fig.  3U.  Block  diagram  manipulations 


for  obtaining  incut-output 
relations  for  servo  system. 


for  obtaining  disturbance 
responses  of  servo  system. 
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These  block  diagram  transformations  are.  effected  by  making  use 
of  the  fact  that  each  component  block  of  the  system  denotes  a linear  operation 
and  the  fact  that  the  overall  operation  performed  on  sampled  input  signals 
by  a continuou3-signal  filter  followed  by  an  impulse  modulator  is  just  that 
performed  by  .a  certain  discrete  signal  filter.  Referring  t.o  Fig.  33,  we 
see  that  the  step  from  33a  to  jjo  results  from  the  linearity  of  the  modulation 
process.  Because  of  this  linearity  it  is  immaterial  whether  the  fed-back 
signal  is  first  subtracted  from  the  input  signal  and  the  result  sampled  or 
whether  both  signals  are  first  sampled  and  then  subtracted.  Step  33b  to  33c 
rests  merely  on  the  observation  that  similar  operations  (here  indicated 
by  the  composite  transfer  function  performed  on  the  same  signal  results 

in  equal  signals.  Step  33c  to  33d  is  the  reduction  of  the  cascade 
filter-modulator  system  to  an  equivalent  discrete  signal  filter.  N eta  that 
at  point  A the  signals  tire  entirely  of  a sampled  nature,  since  they  are 
delivered  from  impulse  modulators.  Thus  we  have  the  branch  with 

discrete  input  signals  which  reduces  to  the  equivalent  filter  with  transfer 
function  or  ' as  described  in  Sec.  1.  of  this  paper. 

Finally  to  go  from  33d  to  33e  we  use  the  well-known  reduction  of  a feedback 
loop  with  open-loop  transfer  function  (ILjKp)  and  feed-forward  transfer 
function  of  1 . 

The  transformations  involved  in  fig , 3U  proceed  similarly.  All 
steps  depend  on  linearity;  in  addition  stap  3hc  to  3Ud  make"  use  of  the 
overall  transformation  given  in  Fig.  33- 

If  we  examine  the  final  result  of  these  manipulations  we  see  that 
we  are  left  with  series-parallel  combinations  of  the  three  fund*anental 
8ampled-data  system  components,  modulators,  discrete-signal  filters,  and 
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continuous-signal  filters.  The  equivalent  systems  obtained  are  the  some 
as  these-  derived  by  use  of  the  more  general  reduction  procedure.  It  is  now 
possible  to  solve  for  the  response  of  these  networks  in  v.  stepwise  manner 
by  means  of  methods  of  Sec.  1.  The  implicit  relationships  in  the  original 
feedback  problem  have  been  reduced  to  expliett.  ones, 

3.  jfflSSTflTPSTQH  0?  DZSIRABIs  aSKV 0 SYSTSH  RESPONSE  CnARACTfeRiSTTSS 

The  earlier  sections  of  this  paper  have  dealt  with  analysis  of 
sampled-data  systems.  Having  thus  completed  analysis,  we  now  turn  to  design. 
Design  is  basically  different  from  analysis.  Analysis  problems  usually 
have  unique  solutions.  Often  there  exist  no  ’unique  design  problems  for 
a practical  situation,  and  even  when  a unique  design  problem  exists  it  seldom 
has  a unique  solution.  At  the  beginning  of  a design  problem  one  has  certain 
design  objectives  on  one  hand  and  certain  physical  limitations  on  the  other. 
At  first  neither  is  clearly  defined.  The  process  of  design  is  a cut— end-try 
process  of  realizing  the-  design  objectives  within  the  framework  of  the 
physical  limitations.  Usually  design  objectives  are  limited  and  made  more 
specific  because  ox  knowledge  of  physical  limitations  and  because  of 
the  designer's  limiting  himself  to  a specific  class  of  situations.  The 
essence  of  evolving  a design  technique  is  to  restrict  and  clarify  both 
objectives  and  realizability  conditions  so  that  the  cut-and-try  process 
can  be  a manageable  one.  Once  evolved,  a design  process  has  three 
component  parts : (l)  formal  statement  of  objectives,  (?)  statement  of 

conditions  of  physical  realizability  conditions,  (3)  an  approximation 
problem  where  objectives  are  realised  within  the  framework  of  realizability 
conditions.  The  last  section  of  this  paper  is  devoted  to  the  study  of  design 
techniques  for  sampled-data  systems.  The  first  step  is  to  present  the 
procedure  for  formal  stater<»ent  of  design  objectives. 
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A sampled-data  servo  is  usually  designed  to  have  specific 
response  characteristics.  The  design  objectives  for  eampled-data  servos 
can  accordingly  be  3tated  in  terms  of  acceptable  response  characteristics. 

In  short,  the  responses  of  a sampled-data  servo  will  probably  be  acceptable 
ifi  (l)  its  transients  are  short  and  adequately  damped,  (2)  disturbances 
in  the  feed-forward  section  are  adequately  suppressed,  O)  errors  in  following 
smooth  input  signals  are  small,  and  (U)  sampling  ripple  in  the  output  is 
small.  It  turns  out  that  the  nature  of  the  transients  can  be  adequately 
characterized  by  locating  the  poles  of  the  overall  system  transfer  function . 

The  quality  with  which  a single -loop  servo  masts  the  next  two  conditions 
can  be  approximately  determined  by  use  of  error  coefficients.  Transient 
studies  made  in  Sec.  1 provide  all  the  tools  necessary  to  analyte  transient 
responses  and  ripple  of  a system  adequately.  Accordingly,  the  response 
characterizations  will  be  completed  by  a study  of  error  coefficients. 

3,1  Error  Coefficients  for  Sampled-Data  Servo3 

Usually  the  signals  handled  by  servo  systems  are  smooth,  with 
only  infrequent  discontinuities.  This  smoothness  of  the  input  signal 
allows  them  to  be  resolved  into  Taylor's  series  much  more  conveniently 
than  they  can  be  resolved  into  Fourier  series  or  into  sets  of  impulses. 

Therefore  it  is  often  better  to  characterize  a servo  system  by  its  error 
response  to  a constant,  a linear  input,  a parabolic  input,  a cubic  input, 
and  so  on,  than  it  is  to  describe  it  in  the  frequency  domain  or  by  its 
impulse  response  alone.  It  should  be  understood  at  the  outset  that  tilts 
error-coefficient  technique  is  merely  another  form  of  superposition  which 
is  particularly  useful  with  servos  because  their  inputs  are  usually  very 
conveniently  resolved  into  power  series.  Furthermore,  servos  arc  most  often 
follow-up  devices,  and  the  objective  is  that  they  shall  follow  with  small  error. 
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The  response  accordingly  can  be  judged  by  knowing  the  error  rather  than 
hv  knowing  the  output.  The'  'successive  error  coefficients  of  a servo 
system  are  the  steady-state  values  of  th»  error  response  at  t = 0 to  a 
constant,  a straight  line  passing  through  the  origin  with  unit  slope,  a 
parabola  passing  through  the  origin  with  unit  second  derivative,  a cubic 
curve  passing  through  the  origin  with  unit  third  derivative,  and  so  on. 

The  great  majority  of  servo  systems  can  bs  described  adequately  in  terms 
of  their  error  coefficients  and  in  terms  of  the  duration  of  their  transients. 

As  with  continuous-data  servos,  sampled-deta  servos  can  also 
be  described  very  well  in  terms  of  error  coefficients.  The  size  of  the 
error  ccsfficients  determines  how  well  a given  servo  will  perform  in 
many  applications.  Since  the  error  in  a sampled-data  servo  is  sampled 
and  since  all  inputs  to  it  are  sampled,  the  error  coefficients  will  be 
found  in  terms  of  sampled  inputs  and  sampled  error.  The  block  diagram  which 
relates  the  input  to  the  error  may  be  obtained  by  the  techniques  of  Sec.  2. 


Fig.  33>.  Equivalent  cascaded  system  of  a one-loop  sampled-data  servo. 

* * 

Tne  relation  between  »Ol  and  can  be  obtained  by  the  transfer  function, 

n a 


1+ 


ZZ* 


, of  the  discrete  filter.  Knowing  the  transfer  function  as  a 


oamee. .Nichols,  and  Tblllips,  Theory  of  Servomechanisms.  McGraw-Hill  Book  Co., 
New  York,  19u7>  Vol. .25-  H.x.T.  Rad.  Lab.  Series. 
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rational  function  of  £ one  can  find  the  unit  sample  response  by 

expanding  the  transfer  function  in  a power  series  of  £ ~ . 

me  positional  error  coefficient  is  the  error  arising  at  t « 0 
from  an  infinite  string  of  unit  sairples  which  have  been  applied  since 
t « - on  . Tnis  error  can  be  evaluated  by  applying  tbs  superposition 
summation  of  En.  7 as  follows:  First  expand  the  discrete-filter  transfer 

function  into  a power  series, 

1 _ -*>T 
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men  apply  Eq.  7 for  the  case  where  n«U  and  Eq.  U»  results* 
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UQ(t)  is  a unit  impulse  at  the  origin.  Thus  Eq,  the  positional  errcr 
coefficient  is  expressed  by  Eq.  !*£  from  the  area  of  the  error  impulse  at 

t-C. 


- area  of  (o) 
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- j.  or 


s =*  u. 

To  find  the  velocity  error  coefficient,  convolve  the  unit  sample 
response  with  an  input  which  is  made  up  cf  the  samples  of  a signal  of 
unit  slope.  For  this  input,  (nT)  * an  impulse  with  area  - nT.  In  Eq.  1*6 
Eq.  7 is  applied  to  the  situation  involving  the  new  type  of  input. 

yt), 

(I46) 

where  UQ(t)  is  a unit  impulse  at  the  origin.  Thus  £^,  the  velocity  error 
coefficient,  is  expressed  by  Eq,  li6  from  the  area  of  the  error  impulse  at 
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The  acceleration  error  coefficient  can  be  obtained  by  applying 
Eq.  7 to  the  case  wherein  S.  (nT)  - an  impulse  with  the  area  ( tiT ) . 


9a(0) 


* m jr  (~n\  * i-2T)^h^  ♦ (-kT)2^  ♦ 


U (t). 


(uo) 


whei-e  UQ(t)  is  a unit  impulse  at  the  origin.  Thus  E^,  the  acceleration 
error  coefficient,  is  expressed  by  Eq.  U8  from  the  area  of  the  error 
impulse  at  t « 0. 


“ area  of  9^ 


(°;  |j  f-w)\  ■ It  (— T£-;i 

k-0  " " T 


s-0 
v 

(1*9) 

The  procedure  indicated  above  can  be  generalized  to  apply  to  all 
error  coefficients  as  in  Eq.  50. 

: / , \ I 

(60) 


R . 1 dk  / 1 \ 

^ *I_  d?  il  ♦ JUkJ 


Note  that  if 


JL 

JT' 


s®0 

<5'sT-l. 

is  expanded  in  a Taylor's  series  about  the 


i<4  sr* 

origin  in  the  s-plane,  the  successive  constants  in  the  Taylor's  expansion 


correspond  to  the  error  coefficient*.  If 


1-42:1 


is  given  in  terms  of  a 


rational  function  of  ^ , it  is  usually  easy  to  make  a Taylor's  series 

of  — ^ — in  terms  of  { ”8^  about  £’"’Sa  » 1 and  then  to  expand  this  series 


l4HK 


✓ .8x 

into  a Taylor's  series  in  s by  expanding  £ - 1 in  2 power  series  in 


1 ♦ f S 


1 ♦ E 


a/  f £ _Si  - \ - -J  t f-  sr  ..  \ 

5 “ o ~1(  ' ’X  i * ’2v(  "1)  * 5 

r _ — 3'T  ..s  _ -sT  ♦ (-sT)2  ^ (-sT)3  . 

(e  -1)  - -jj * + 


(6D 

(52) 
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Combining  Eqa  - ?>1  and  52  yields  2a.  53  < 


r..* 


®f  . ♦ i ♦ -W 


(-sT)\  (-eT)3  A 


1 ♦ K 


-H 

~ / (-sT).  C-sT)^  (-aT)-'  ^ 

^2  ”71“  “TP  """ 


From  Eq.  53  it  can  be  seen  that  Eo*°^o(-T)0,  E^» 

Eo-  &S  WJ(-T)2,  etc.  Usually  the  positional  error  coefficient  is 

L.  C.  C X 

zero,  and  often  the  velocity  error  coefficient  ie  also  zero-  Since  ve 
are  most  interested  in  the  first  few  error  coefficients  and  particularly 
in  the  first  non-zero  one,  the  Taylor  series  in  (6~81-!/,  Eq.  5l,  is  usually 
the  most  convenient  relationship  to  use.  Suppose  that  — — — — has  a 

1 ♦ r" 

kth  order  were  at  “ _ST-  1 '”sT'.  i~.~. 


1 ♦ K 


“J.  / f V ^ ) * All 


of  the  first  (k-l)  error  coefficients  will  be  zero  snd  will  bs  given 
by  (-T)kF(  £ 8l)  | . If'F(£”s^)  is  a rational  function  of  ^ , it  can 


t -1 


be  evaluated  by  taking  t.ho  rctio  of  the  product  of  all  the  vectors  whose 
tills  rest  on  zeros  of  F and  whose  tips  rest  on  the  >1  point  to  the  product 
of  all  the  vectors  whose  tails  lie  on  th®  poles  of  ? ana  whose  tips  rest 


on  the  ♦!  point  <.  Thu 3 a zero-velocity  error  nerve  must  have  a 


»5T 

function  having  a double  order  zero  at  r =1.  It  is  interesting  to  note 

in  passing  that  the  behavior  of  a sampled-data  servo  corresponds  to  the 

-sT 

behavior  of  its  transfer  functions  about  q »1  in  the  same  sense  that  the 
behavior  of  a conventional  servo  corresponds  to  the  behavior  of  its  transfer 
function  about  s-C. 


i 
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Having  found  a procedure  for  evaluating  the  error  coefficients,, 
one  can  apply  them  to  evaluate  the  error  resulting  from  a lead  torque 
disturbance  much  as  is  done  with  continuous  servos.  3y  methods  shown  in 
Sec.  2t  a steady  load  torque  disturbance  can  be  seen  to  cause  the  same 
steady  error  as  a velocity  input,  and  consequently  low  velocity  error 
coefficients  insure  good  suppression  of  load  torque  disturbances. 

Error  coefficients  ana  pole  positions  of  the  system  function  are 
sufficient  to  describe  the.  essential  response  characteristics  of  the  system. 
Error  coefficients,  and  the  nature  of  the  transient  response  can.  be  obtained 
from  a study  of  the  i function.  Accordingly,  from  now  on  the 

response  characteristics  will  be  characterized  by  the  constellation  of  poles 
3 -sT 

and  zeros  of — in  the  C -plane.  It  should  bs  emphasised  that 

s . i-  V '*■ 

A- 

this  study  of  the  sampled  signals  does  not  give  a measure  of  the  ripple 

in  the  continuous  output,  but  that  this  measure  car.  easily  be  obtained 

in  a specific  case  from  an  analysis  such  as  was  giver,  in  osc . 1.  Since 

the  design  objectives  must  be  cut  to  the  bone  to  make  a design  procedure 

workable,  it  seems  proper  to  assume  that  if  the  sampled  output  transients 

are  acceptably  smooth  and  if  the  output  section  is  low'-pass  enough,  the 

continuous-transient  response  will  be  acceptable.  Always  after  a given 

design  is  worked  out  the  ripple  c an  be  calculated  for  specific  cases  * 

3*'?  Restrictions  Placed  or.  Compensation  bj  Fbysiiai  Realizability 
Conditions 

The  previous  section  has  characterized  by  a somewhat  arbitrary 
(though  sensible)  limiting  of  the  desired  response  description.  This  section 
is  characterized  by  » similar  limiting  of  the  types  of  compensation  considered. 


The  procedure  presented  is  applicable  to  a wide  variety  of  problems.  In 
principle  it  may  be  extended  to  fit  a much  wider  range  of  situations.  The 
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point  of  view  and  the  analysis  teehnaiue  are  actually  more  important  than 
the  formal  compensation  procedure  specifically  presented.  The  reader  is 
again  reminded  that  to  be  workable  a design  procedure  must  be  much  less 
than  completely  general.  The  design  procedure  given  here  would  serve  as 
a starting  point  for  the  study  of  a specific  design,  and  s.3  the  design 
progressed  a more  elaborate  s«-t  of  procedures  could  be  worked  out.  As  with 
continuous  servos,  practical  considerations  might  make  the  specific 
configuration  presented  here  unworkable.  Many  of  the  conditions  derive' 
here  would  carry  over  to  other  situations,  so  the  comments  below  are  move 
general  than  the  situation  they  are  applied  to.  Figure  36  shows  the 
m compensated  system  and  the  type  of  compensation  proposed.  The  system 


(a)  Uncompensated  system  (b)  Compensated  system 

Fii'.  36.  A sisiple  servo  without  and  with  compensation  in  the  forward  loop. 

has  unity  feedback,  and  the  compensation  filter  can  be  an  entirely  discrete 
filter,  an  entirely  continuous  filter,  or  a combination  of  the  two.  The 

overall  input- to -error  transfer  function,  — -i- , is  the  one  used.  It 

is  modified  by  compensation  networks  to  change  it  to  more  acceptable  values. 

The  oitput  element  has  certain  characteristics  which  can  be  changed  only 
at  considerable  expense.  Accordingly  they  are  assumed  fixed.  The  compensation 
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technique  is  to  select  a desirable  — ■ — =, function  and  the!1.  w ork 

backward  to  the  compensating  filters.  First  physical  roaiizability  conditions 
should  be  placed  on  the  filter  transfer  functions.  From  these  filter 
realizability  conditions,  conditions  on  derived. 

3.21  Physical  Realizability  Conditions  on  Transfer  Functions  of 
niscrsts  r.nd~Sontinuous  Filters 

The  realizability  conditions  on  transfer  functions  of  a Bingle 

filter  are  fairly  simple  to  state.  For  a discrete-signal  filter,  the  transfo* 

function  must  be  a rational  function  in  £ *<>  The  fundamental  realizability 

condition  for  any  filter  is  that  the  response  cannot  occur  before  the 

input.  This  condition  dictates. that  the  t roister  function  have  no  poles 

— gT 

at  the  origin  in  the  "'-plane:  it  is  t.he  only  realizability  condition 

on  discrete  filters  which  is  fundamental.  Usually  a discrete  filter  would  • 

be  expected  to  be  stable,  though  unstable  filters  can  be  realized.  A 

stable  discrete  filter  must  have  transfer-function  poles  located  only 

••sT 

outside  the  unit  circle  centered  at  the  origin  of  the  £ -plane.  The 
fundamental  physical  realizability  condition  for  corMrvous  filters  is 
that  +-he  response  cannot  occur  before  the  excitation,  and  this  requirement 
is  translated  into  the  frequency  domain  by  the  Paley-Wiener  criterion. 

Actually  we  will  deal  with  transfer  functions  which  are  rational  functions 
of  s which  all  satisfy  the  Paley-Wiener  criterion.  For  practical  reasons, 
such  as  parasitic  capacitance  in  electric  filters,  it  is  not  feasible  to 
have  transfer  functions  having  more  zeros  than  poles  in  the  finite  s-piana. 
Accordingly,  usable  continuous-filter  transfer  functions  will  be  limited 
to  those  rational  functions  of  s having  no  mere  poles  than  zeros  in  the 
finite  s-plane. 
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Having  tie  conditions  of  physical  realizability  compactly  stated 


for  component  filters  which  will  be  cascaded  together  to  make  up  the 
feed- forward  section  of  a sampler-data  servo,  the  designer  must  now  tun 
to  the  problem  of  assossing  what  limitations  are  placed  on  the  overall 


1 v — » 


Trans f : 


function  by  tne  fixed  part  of  the  system. 

3.22  Physical  Reaiizabiiitv  Limitations  Placed  on  Feed-forward 


by  Limitations  of  Fixed  Member 


tr lotion  to  Discrete  Filter  Compensation 


Given  a fixed  output  member  with  a transfer  function  K.(s), 

iJsaa  r 

one  can  calculate  its  sampled  equivalent  as  K^(s»jnXj).  If  the 


— — sT 

discrete  filter  to  be  used  has  a transfer  function  D( £ ),  then  the 

overall  sampled  transfer  function  for  the  feed-forward  section  is 
n “^)  * F_(s'»*r»n).  The  physical  realizability  condition  of 

n'-Oo 

the  whole  sampled  feed-forward  furiction  is  then  restricted  only  to  guarantee 
— bT 

that  D( £ ) be  realizable  or  that  it  have  no  poles  at  the  origin.  In 

short,  the  overall  sampled  foed-forward  transfer  function  will  be  realizable 
"if  It  has  at  least  as  many  zeros  at  the  origin  of  the  £ -plane  33 
K^Cs+JnTl)  does. 


n*—  00 


3.23  Physical  Realizability  limitations  Placed  on  Feed-forward 
Transfer  Functions,  by  Limitations  on  Fixed  Member 

and  Restriction  to  Continuous  Filter  Compensation 


If  the  designer  were  given  a fixed  continuous-signal  output 
member  with  a transfer  function  K (s)  and  if  he  had  selected  a sampled 

1 » 1 00 

overall  feed-forward  transfer  function  X (s)  . K(s+jnXl)  he  would 


find  the  compensating  network  as  follows: 


n™-»o 
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(l)  From  ^ \ K's+jnn)  he  would  find  a suitable  K(s).  (One 
n"-  ®o 

K (s)  always  may  correspond  to  any  one  of  several  K(s) 

£. 

functions,) 

(?)  The  transfer  function  of  the  compensating  network  would  be 
K/Kf. 

at 

The  problem  of  this  part  is  to  find  how  K~(s)  is  restricted  by  characteristics 


of  K*(s)  and  the  fact  that  ~ 
x K. 


must  be  physically  realizable.  As  the  next. 


few  paragraphs  will  show  the  only  conditions  on  K (s)  such  that  K/E^  be 
physically  realizable  are  these* 

Cl)  If  a.i's)  has  an  excess  of  poles  over  zeros  of  two  or  more 

“ rsT 

K f s)  must  havs  a single  zero  at  the  origin,  in  the  £ — plane* 

(2)  K (s)  must  have  one  less  zero  than  pole  in  the  £ -plane 

so  that  it  can  be  expanded  in  a partial  fraction  expansion 

with  no  constant  term* 

Condition  (l)  results  from  the  fact  that  « — can  have  no  more 

Rf 

zeros  than  poles  in  tha  finite  s-plane  and  accordingly  K must  have  the 

same  excess  in  numbers  of  poles  over  zeros  in  the  finite  s-plane  that  K„ 

* l^s^1 

has.  For  5 to  have  two  more  poles  than  zeros  requires  that  K « ~ y K(s*jn-0) 


-sT 


n»-  c*o 


must  have  a simple  zero  at  the  origin  in  the  ^ -plane,  but  if 


IT.  v 

K (s)  has  a single  zero  at  the  origin,  an  equivalent  unsampled  K can  be 
found  having  any  finite  number  more  poles  than  zeros.  One  possible 
procedure  will  be  outlined  below. 

1.  Resolve  into  a partial  fraction  expansion.  If  all 

' ' " ' ' ' ‘ ' .*  1' 


The  fact  that  a given  sampled  transfer  function  K"-  ^ \ E(s*tjn.Q}  does 

n~~ocj 

not  lead  uniquely  to  a K(s)  is  exactly  equivalent  to  the  fart  that  several 
different  time  functions  can  *11  lead  to  the  same  set  of  sar.pl so  sc  long 
ns  they  are  alike  only  at  the  sampling  points.  Two  filters  having  different 
E(s)  functions  will  have  the  same  K (s)  functions  if  their  impulse  responses 
arc  alitce  at  the  sampling  points  but  different  between  sampling  points. 
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•poles  are  not  simple-  treat  multiple-order  poLes  »3  a group 
of  simple  poles  with  small  saporatien*  The  form  of  the 
expansion  is  given  in  Eq. 


Kt/  \ ! 

is;  x i 


r^. 


- ±<*. 


rr~ - — na — ~ X i 

1_  £-'s-Va  2 1 
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- i I ♦ . 

2 ?.  j 

a5U) 


2.  To  each  term  in  the  partial  fraction  expansion  c&ssoci  fltifi 


an  unsampled  filter  transfer  function. 


-i s-a  )T  ~ * 

i-  v v - 


p. 


(o) 


(-1) 


ft  (1) 


B-I, 


s-a.tjJT. 


♦ . 


*i 
(-k) 


(55) 


The  only  constraint  anorig  the  p.,  s is  that  P4 

(k)  • j 

fh  be  conjugates  and  that  the  sum  of  all  p^  s adds  up  to<»<  j.  » The 

t \ 

simplest  filter  results  from  using  the  low-pass  term  PiVO/  only  for 

, s-a_. 

each  term  ir.  the  partial-fraction  expansion  of  x 

1 T"  £ 

If  Y~y  K has  a zero  at  the  origin;  then  y i“^*  result 

l5T  n 

can  be  seen  by  adding  together  the  terms  in  F.q.  5U-  If  'y  ^-0,  then 

i-1 

n ra (k)  n (o) 

\ \ p *0.  If  only  the  low-pass  equivalent  is  used,  then  \ P ..  * 

i-1  0 ' i-1 

or  the  continuous  filter  has  the  sum  of  all  of  its  residues  equal  t-o  zero 

whence  it  must  have  two  more  poles  than  zeros.  The  additional  freedom 

in  selecting  the  p's  allows  one  to  make  sl(s) — 0 as  ~2  as  s ^ oo 

s 

so  that  K can  be  made  to  have  the  denominator  thr««  degrees  higher  than 
the  numerator.  By  continuing  to  use  the  additional  freedom  in  selecting 
P's,  one  can  make  the  degree  of  the  denominator  of  an  equivalent  I.  be 
as  many  degrees  higher  than  the  numerator  as  he  desires.  The  only  cost 
is  a more  complicated  K(s)  function.  The  procedure  to  select  an  equivalent 
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K with  a denominator  m degrees  higher  than  the  numerator  is  to  select 

the  P's  such  that  the  sum  of  ?11  p.'s  add  up  to  and  further,  that 

the  p!s  have  such  values  that:  the  sum  of  the  residues  ox  csX  equal  zero, 

2 

the  sum  of  the  residues  of  a K equal  zero,  and.  so  on  including  the  sum 

in*-*** 

of  the  residues  of  s equal  zero.  *3  the  degree  of  the  denominator 
exceeds  the  degree  of  the  numerator  by  a greater  and  greater  amount,  the 
number  of  simultaneous  algebraic  equations  the  P's  must  satisfy  increases. 

T).ie  fact  that  & given  If  (s)  having  a zero  at  the  origin  in  the 
(f  -plane  can  be  realized  by  a I(s)  having  a much  higher  degree  denuminaur 
than  numerator  is  physically  reasonabl-s.  If  a E*(s)  function  has  a single 
zero  at  the  origin,  the  initial  sample  of  the  unit-sample  response  will 
be  zero.  The  fact  that  the  degree  of  the  denominator  of  f(s)  is  much  higher 
than  the  numerator  requires  that  many  initial  derivatives  of  the  impulse 


■*'  S.  J OT  T r>« 

y.  -*  W V/J.  VtiC. 


he  continuous  filter  be  zero  but  does  not  limit  what  values 


it  can  have  at  T seconds  after  the  impulse  was  applied.  Accordingly 

g. 

one  would  expect  the  requirement  that  S ( w)  huu  a oCrO  at  bli6  origin  if 
the  degree  of  K's  denominator  exceeds  the  degree  of  its  numerator  by  two 
or  more  and  is  not  surprised  that  there  are  no  other  restrictions  an  the 
sampled  transfer-  function. 

The  restrictions  placed  on  by  Uotj  of  discrete-filter 

compensation  and  by  continuous-filter  compensation  ha  Vc  Sxiua  laritie3.  If 
K^Cs)  has  two  more  poles  than  zero^  both  types  of  compensation  must  lead 
to  iT^K'a  which  have  a single  zero  at  the  origin  in  the  ~ -plane 


Given  1C  (s)  one  can  employ  st^ll  furtner  freedom  in  going  to  a K(s)  function. 
Given  any  K(s)  which  yields  1C  (s),  one  can  add  to  it  any  further  transfer 
function  whose  impulse  response  goes  through  zero  at  all  the  sampling 
points.  One  such  zero-sampie  function  is  P where  p and 

(s+c4)^+  /lf/U 

are  any  two  real  numbers.  See  Appendix  B for  further  details. 
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ieg~x-uless  of  vhetner  discrete  or  continuous-filter  compensation  are  used. 

For  continuous-filter  compensation  the  number  of  zeros  of  ^ K must  be 

one  less  than  the  number  of  its  poles  (Condition  2)  because  there  is  no 

continuous-signal  lumped-parameter  filter  realization  for  a K (s)  function 

which  is  a polynomial  in  £ "3l,  (It  could  be  realized  by  delay  lines, 

however.)  "his  restriction  does  not  apply  to  discrete  filters.  For  discrete 

_ -aT 

filter  compensation  the  number  of  zeros  at  the  origin  in  r~  -plane  of 

1 T 

D - ^ K^.  must  be  no  less  than  the  number  of  zeros  of  This 

restriction  does  not  apply  to  continuous  filters. 

1 < 

In  any  case,  the  compensation  is  simpler  if  the  poles  of 
are  included  among  the  poles  of  than  if  is  selected  without 

attention  tc  the  poles  of  use  of  a combination  cl  discrete  and 

continuous  filters  i3  often  useful  because  it  allows  one  additional 
freedom,  in  getting  K(s)  which  may  then  be  used  for  ripple  suppression. 

The  use  of  a combination  of  discrete  and  continuous  filters  results  in 
more  relaxed  realizability  conditions  than  wou ) J be  required  for  either 
discrete  or  continuous  filters  alone.  The  only  condition  of  realizability 
on  fO  for  compensation  by  a combination  of  discrete  and  continuous 

1 V *' 

i liters  is  that  if  K,  has  two  more  Doles  than  zeros  —> K must  have  a 

I 1 

/ ■g1]’ 

zero  at  the  origin  in  the  c -plane.  Aside  from  this  restriction, 

-sT 

any  rational  function  of  f is  realizable. 

Note  that  if  K (s)  Includes  the  transfer  function  of  a clamper  then  the 
degree  of  the  numerator  of  a (s)  can  be  equal  to  the  degree  of  the  denominator 

Can  ku)  - TjgSsl].  K*.  (1_e-»V| 

Given  K (s),  to  find  — - one  would  resolve  — jtT.  into  e part.;.!  fraction 


expansion.  Such  an  expansion  would  have  no  constant  term  if  the  decree 

of  the  numerator  of  — hit,  vc?  1 less  than  the  degree  of  the  denominator. 

_ ~ ^ 
mis  requirement  demands  that  the  degree  of  the  numerator  of  K (s)  be 

no  greater  than  the  degree  of  the  denominator. 
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From  the  description  of  realizability  conditions  on  the  sampled 
feed-forward  transfer  function  one  can  deduce  conditions  of  physical 
realizability  on  the  function  — sp 


1*  so 


3.2U  Restrictions  on 


1 

"T 


IT~K 

~ T, — 


Function  from  Realize!  iiity  Conditions 


on  the  Compensating  Network . 

Given  the  realizability  conditions  on. the  feed-forward  transfer 
function,  the  conditions  on  — r can  be  derived. 

Hn* 

Lett 


K*(s)  - K.Ca+Jnfl)  * ■ - m‘ ' . 


(56) 


n--< 


i + iTr. 


F*Q 


i \ 

i .Of  ) 


Given 


-sT 

as  a rational  function  of  ^ *,  one  immediately 


1 Nr-*-  ~ 

1+T2-* 

has  Q from  the  numerator  polynomial.  Knowing  Q and  (P*Q)  one  finds  IJ  by 


finding  the  difference.  P will  }i»v»  a zero  at  the  origin  of 


■ M 


-sT 


(58) 


.-sT 


The  design  in  ten?; s of  -p^Q—  is  particularly  easy  because  setting  the  zeros 
of  the  fraction  at  the  polos  of  the  sampled  fixed  function  leads  to  simple 
compensation  networks. 


! 
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3-3  The  Design  of  oanpled-Data  Systems  by  Approximating  to  Desigr 
Ob.i ectives  with  Realizaole  Compensation 

The  design  procedure  sketched  below  is  somewhat  aribtrary  and 
my  be  edified  to  adept  it  to  specific  situations.  The  steps  of  the  design 
procedure  are  these: 

1.  Analyze  the  fixed  parts  of  the  system.  Find  the  position 

1 1 “ST 

of  the  poles  of  -^ZKf  in  the  (f  -plane  end  make  sure 

that  the  sampling  frequency,  the  bandwidth  of  the  servomotor 

elements,  and  the  bandwidth  of  expected  signals  are  consistent. 

1 


2.  Sketch  a map  of  desired  poles  and  zeros  of 


i*$Ek 


in 


^Ciiwo  HUSLR.t3 


the  £ “ -plane.  Evaluate  the  error  cceffici 
sure  if  Q is  to  have  a zero  at  the  origin  that  the  r>roduct 
of  pole  distances  from  the  origin  equals  in  both  magnitude 
and  angle  the  product  of  zero  distances.  Whether  the  system 
transients  are  proper  is  dependent  largely  upon  the  location 

i 

of  the  dominant  poles  of 


3.  Having  selected 


1 p 

, solve  f or -j  and  then  select  a 


l.f£K 


D and  a K which  realizes  ^ and  a«.  the  same  time  provides 
proper  smoothing  nf  sampling  i-ipple.  This  process  is  largely 
cut-and-try  where  one  tries  the  simplest  realizations  first 
and  proceeds  to  more  complicated  realizations  as  necessity 
demands . 

Since  the  body  of  techniques  for  ssmpled-data  servo  system  design 
is  quite  extensive,  it  is  not  feasible  to  catalogue  ail  of  them,  but  rather 
to  indicate  that  most  usual  techniques  for  ccntiiiuous-data  servo  system 
design  have  counterparts  in  the  sarapled-data  ease  and  then  give  a few 
i] lustraid  examples . 
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3.31  Design  of  a 5°mpIcd=Dsta  Distance  indicator 
A distance  measuring  device  is  to  be  made  to  indicate  both 
distance  and  velocity  as  continuous  quantities.  The  distance  of  the 
object  is  provided  by  a puls exposition  modulated  signal  with  a repetition 
frequency  of  one  pulse  per  second.  The  mood. mum  pulse  delay  is  u srcall 
fraction  of  a second.  The  block  diagram  of  the  device  is  shown  in  Fig*  3? 
below.  The  practical  problems  associated  with  the  design  of  the  system 


Voltage  proportional  to 
velocity 


P-I'-M 


Dr  ] fty 

direct  voltage 

Compensating 

_J 

Integrator 

Detector 

to  delay  beiweei 

1 Network 

- 

and  fed-back  pulse 


Pules -Position 
Modulator 


Voltage  proportional 
to  the 


(a)  Block  diagram,  of  equipaeni 


Fig.  37-  The  distance  measuring  device  and  its  equivalent  block  diagram* 
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(i)  Tins  delay  detector  has  a very  limited  dyncnic  range  and  is 


linear  only  aro^n^  i+y  null  point.  (2)  There  i*  some  noise  in  tho  actual 
distance  indicating  delay  so  that  while  the  pure  inputs  to  the  device 
would  correspond  to  sets  of  rampe  actually  there  are  superposed  noise 
signals  having  bandwidths  of  cycle  par  second  cr  less.  (3)  The  output 
should  follow  the  actual  signal  as  closely  as  possible,  (i*)  The  gain  of 
the  null  detector  has  some  drift  associated  with  it.  These  gain  changes 
should  not  adversely  affect  the  system  performance. 

The  first  step  in  the  design  procedure  is  to  evaluate  and 


K_(s)  - i- 


K - ( a ) 


a-  r8T)  \ 


f (a4.1n«f  L 

n«—  o>© 


(i-  rsT>?  raT ..  t f 

zw* 

(l-  e SV 


i-  f 

(56) 


.-ST  * 


* 

Sines  hae  « icry  •*!  the  origin  in  the  ^ -plane,  the  overall 

feed-forward  function  must  have  this  same  «rc.  This  is  the  only  limitation 

1 


on  R . The  selection  of  a 


function  can  now  be  made.  The  selection 


1+K 


of  the 


1*K 


function  if  by  no  means  a unique.  The  more  data  the  engineer 


has  on  the  specific  problem,  the  more  intelligently  he  can  use  the  freedom 

* 

in  selecting  K be  has  gained  from  knowledge  of  general  realizahi'lity 

is  the  best  function  to  design 


* 1 

conditions  on  K . In  many  cases 


1+R* 

in  terms  of  but  it  is  not  the  only  one.  In  this  exanple  we  will  make  a 

fairly  arbitrary  selection  of  K and  go  through  the  design  before  looking 

fer  the  alternative  possibilities.  In  order  to  guarantee  a zero  velocity 
1 - --*T 


# must  have  a double  order  zero  at  the  ♦!  point  in  the  ' -plane 

i+r 

for  cur*  type  of  condensation.  To  make  the  transient  short  move  all  the 

1 

l*t 


error, 

1*K" 

4- 

''Jb'V 

poles  of  tho 


function  to  irfinitj*  A3  a first  trial,  assume 
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n« 

-OX- 


1*K  

Since  the  denominator  of 


r,  /-sTv2  « , -sT,  1 -six 

-•  ' — * ‘ — * The  resulting  value  of  K is  ^ A 2 ^ 

- L 1 ~^Tsr7 — 


and  its  numerator  have  the  same  value  at 


1+S 


-sT  * 

T “0.  the  numerator  of  K is  ernarantoed  to  have  a root  at  the  origin 

44 

The  cons  tell  at?.on  of  poles  and  zeros  of  IT  is  plotted  in  Fig.  3d» 


-sT 


-1 


r ae 


X - 2 /3 


^ I « 2 


locus  of  -1 
value?  of  the 

_* 

a junction 


-{plane 


J-K  - : 


I A 
i 

5 


J meo 


? - 22/3 


1 

Fig.  3S.  Poles  and  zeros  of  K in  the  £ -plane  and  locue  of  poles 

n* 

cf  the  system  function  *o  when  the  gain  is  changed. 

°1# 

Note  if  the  gain  of  T is  1 the  -l’s  of  I*  occur  at  £“sl»  li  3:  if  E-2. 

& 2 * 

tli*  -l1?  of  X occur  ot  infinity*  If  X-  2 /3*  the  -1*3  of  K occur  at 

x — sT 

C “ -1  and  *3.  This  system  has  a very  good  transient  response;  a zero 
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velocity  error  coefficient  and  the  system  stability  has  fairly  high  eenattivity 
to  changes  in  loop  gain.  Note  that  the  steedy-state  response  to  steps  or 
rasps  still  has  a z-fro  error  far  wide  changes  in  E.  Whether  one  could 
tolerate  the  high  sensitivity  of  the  system  stability  to  changes  in  loop 
gain  would  be  dependent  upon  tne  variability  in  gain  of  the  error  detector. 

file  iiSjli/  u 1/1 ^ . r whs  design  is  to  calculate  the  transfer  function 

of  the  filter  needed  to  obtain  the  compensation  chosen. 

K*  - 12 -1  1 ? 1 . (Si) 

a-  rsT)2 

, . ' 1-  ..  . , 

Pince  tne  clamper  ties  a uransier  luncoxon  * unen  ons  rjiamper 

8 

discrete  filter  part,  (l-  ^ ) must  be  removed  and  the  continuous  filter 

must  give  a K,  as  indicated  by  Eq.  58. 

«*.  *rsT<i-ir\  (53> 


. , -sT\j 

(1-  t ' 


mt  jc x xi x » " U . ..  .A  a*l  X U - X 11. 

ills  XdLUu  uiido  ii cio  ri  uaxi\*  puxw  tx o uio<&iio  iuoo  ui*S  CuaT/iTiuuuc 

filter  iuUst  have  a third  order  pole  at  the  origin.  From  table  A2  in  the 


appendix  note  that  £ / - , „ __  , . - 

A ^(^nn)^  2 a-  rBT)> 

s 

is  to  subtract  from  a multiple  of  the  function  above  such  that  the 

. -«T 

remainder  has  no  third-order  pole  at  ^ **1. 


*2  c-zT , 


• The  probless 


-*(!♦  6”?t)  - 2 f sTci-  \ e*‘sT) 


AT  - 1, 

T 

A - -r 
T 


r8T-l 


I ro*r  ?4  j f 

(l-  £ _sT)3  Cl- 
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,rhe  continuous  filter  must  nave  the  transfer  function  given  by  3q*  6.1 . 


K_  - 


(6i) 


The  compensating  network  transfer  function  is  qiven  by  dividing  K by 


K comp.  “ 


, ,3T 

vsr— 

£ 8 


(62) 


Fig.  3 ?.  The  proposed  system  with  the  conpensatlng  network. 

Having  selected  a system  function  — » — -r — and  having  seen  how 

1*K 

it  can  be  realized,  let  us  now  investigate  further  what  the  response 
characteristics  of  the  system  are.  It  is  easy  to  get  the  sampled  impulse 
response  * Call  the  outpxit  of  the  inpulse  modulator  ft  (t)  (Pig*  38). 

St 


1*K 


i - 2 -T.  e -2sT 


(63) 


r\-< 


If  y^(t)  is  a single  impulse,  then  the  response  can  be  seen  as  in  fig* 

£. 

For  S4(t),  a sartp?Led  step,  the  response  is  seen  to  be  a step  starting 

J. 

wm  rt  + Vi  A (KMMnl  A A<*  O «+  +■  = *£  «*»/<  «A  vVlO  *»P  oft  pr  . 

O O 4 V.  HxOtl  WUHlj/^-w  w*  u M v w /»  vwtwi  <<v  v*  • v*  .*»  Vi 'JA  OCi  * 


for  y.  a unit  impulse 

ft  1 

£ * * 

T for  a unit  impulse 


\ 


£•  & 

eg  for  0.,  a unit  step 


r y* 

f :q 

I T 


f ir  ©_.  a unit  ramp 


J i 


4 

* 


3T 
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Fig.  U0  Sampled  outputs  for  several  inputs. 

Note  that  the  step  response  has  an  overshoot  such  that  the  sum  of 
errors  is  aero.  This  overshooting  characteristic  is  a necessary  part 


ho. 

fit- 


6T  t 


..>t-eV3» 
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'-0>» 


of  the  zero -velocity-error  characteristic.  Two  xcuujj  x-eaponse  is  piouoea 

in  fig.  hO-  The  ramp  response  can  bo  calculated  from  the  step  response 

by  a very  simple  procedure.  Note  that  the  ramp  input  samples  car 

be  obtained  readi3.y  from  the  unit  step  input  s r moles  by  delaying  T.hs-. 

by  1 sampling  interval  and  then  multiplying  by  T and  adding  up  ail 

the  present  and  all  previous  samples.  Thus  the  ramp  at  t - T equals 

T or  the  T times  the  sample  at  t - 0 plus  all  previous  ones  (or  zero). 

The  rarap  sample  at  t - 3T  is  T times  the  sum  of  step  samples  at 

t = 2T,  T,  and  0 or  3Te  Similarly  the  response  to  the  ramp  at  t ■»  3T 

is  found  by  multiplying  by  T the  sum  of  the  step  responses  at  t ■ 2T, 

T,  0.  If  there  is  to  be  no  ramp  error  then  the  sum  of  the  step 

rsspone  samples  must  equal  the  sum  of  the  step  input  samples. 

Having  a picture  of  the  samples  of  the  ramp  response,  next 

we  sneuid  try  to  get  a measure  of  the  actual  continuous  output. 

The  til*  fins f* orm  the  arf-USi.  >fh.9n  ° TSXUp  input  is  Sppll;';d  il* 

given  by  multiplying  the  transform  of  the  sampled  ramp  by  the 

transfer  function  relating  (s)  to  *(s) . 

c X 


^(s) 


-sT 


(l-  e sT)2 


r~i  s-  -sT.2 
O-  c ; x 


Ci-  £“bT)(i*  2jr  •> 

“77 


C6U) 


© (s) 
c 


(l-  6 


-sT-, 


_«T. 


(!♦ 


Ter* 


3T 

2s) 


(65) 
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Since*9^(s)  has  nc  poles  except  at  the  origin,  it  can  readily  be  seer,  that 
the  only  types  of  long-term  response  are  parabolas,  roups  or  steps. 


/ - 

'‘t(s>  ■ e x(r5 

Vis 


1 1 . ^-2s v: 


2/ 


•Tr- 


,66) 


The  output  time  function  can  be  evaluated  by  evsulating  the  inverse  transform 

1 12 

of  the  elements  of  G_  . The  inverse  transform  of  — *•  is  -^p  i for  t greater 

?s 

than  zero.  The  invert*  transform  of  ^ -rj-  - ^ t for  t greater  than  zero* 

6 

The  inverse  transform  of  ”sT  ^ is  "vbiCt-T)2  + 4 (t-T)  for 

Ts3  d s"  2 

t>T.  The  inverse  transform  of  ^”2s'  /™^T*  4 "^y  1 is 

\Ts'  ^ s"  / 

- ^ — (t-2T)  - — Ct-2'i’)  for  t >2T.  At  all  other  values  of  t the  response 

is  zero.  For  t _>?T  tbe  total  response  is 


f(t--T)2  - pt-T)-fjH  ♦ \ jt-T  - t*2TJ  = (t-T)  + 


T.  ■ t. 


The  output  curve  is  plotted  in  Fig.  id.. 
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Having  followed  through  ,v.  trial  design  we  should  new  observe  ir 
retrospect,  whet  was  done.  In  the  first  place  the  process  illustrated  here 
was  not  a.  ‘'one-shot"  synthesis  procedure  to  pet  an  "optimum  system."  The 
process  is  trial-and-errcr  and  what  has  been  done  so  far  is  to  make  one 
trial  design  and  evaluate  it.  The  fact  that  system  characteristics  are 
given  compactly  in  terns  of  the  samplea-irput- 1 o -sasspled-error  transfer 
function  and  that  realizability  conditions  are  easy  to  -state  means  that  the 
first  trial  can  be  made  with  enough  circumspection  to  be  a successful 
on  6 • r or  cm*  6 iwcuiipl.6  the  first  trial  could  be  characterized  as  possessing 


very  good  steady-state  characteristics  and  short,  transients*  The  fact 
that  there  is  high  overshoot  to  the  step  resnense  means  that  the  system 
would  not  be  satisfactory  to  handle  1 ..put - uith  many  discontinuities 
in  them,  Response  to  a ramp  is  certainly  smooth  and  fast.  Fig,  33 
indicates  that  the  system  will  not  be  overly  sensitive  to  gain  change* 
if  discontinuities  in  the  input  occur  often  enough  so  that  the  high 


overshoot  in  the  step  response  is  objectionable,  then  the  — could 

1+K 

be  rc chosen.  The  fact  that  a zen: -velocity-error  system  is  specified 
implies  that  the  step  response  must  overshoot  if  it  starts  off  at  zero. 


but  the  severity  of  the  overshoot  can  be  diminished  by  making  the  cver- 

1 


shoots  Xast*.  !L on Xif 


a-  e’"sT'>2 

■ - « ■ ^ whv  po&rC  u V 6I'31100X» 


i+k“  , i-  rsT/2 


to  a steo  can  be  seen  to  bo  ■?;  but  the  transient  lasts  for  a much  longer 


time.  It  should  be  observed  that  one  has  much  freedom  in  the  design 
procedure  but  to  use  it  intelligently  he  must  have  many  practical  details 
of  the  system  in  mind.  Since  the  problem  chosen  hero  was  merely  used 
as  an  example,  it  will  not  be  carried  further. 
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3.32  Compensation  of  a Sampled-Data  S/stem  with  a Motor  Output 
Consicer  the  system  pictured  in  Fig.  u2.  Suppose  the  sampling 


Fig.  42..  A unity  feedback  3ystem  with  a motor  output. 

interval  is  1 second  and  that  the  motor  time  constant  7~  is  0.?21  second. 

design  a sys terra  vo.th  sere  velocity  error,  short  transients,  and  one  which 

requires  only  a continuous  .filter  for  compensation.  (Actually  one  could 

prescribe  a discrete  compensating  filter  but  it  would  be  more  expensive 

to  realize  than  a passive  electric  filter  would  be.)  The  transfer  function 

of  the  feed-forward  section  is  (l-  f -sT")— *— K fs).  The  (l-  ^-3^^ 

'a^ZVl)  C 

factor  is  essentially  that  of  a discrete  filter,  and  the  fixed  continuous-filter 

1 


transfer  function 


s^T^s*!) 


implies  that  the  overall  continuous  filter 


must  have  uhree  more  poles  than  zeros.  The  pole  at  s*  - a implies  a pole 

* T * * 

of  K at  £ “ U.  From  these  restrictions  on  K one  can  select  reasonable 

1 


1+K 


functions.  As  a start  try 


(l-  6 ~sT)  (l-  esT/u) 


lo  I) 


1*K  1 

_ m 

This  insures  a zero  velocity  error  system  with  a zero  at  £ to  use 

fr 

the  pole  of  K arising  from  the  fixed  part  of  the  system.  The  impulse 

1 


response  of  the  system  is  obtained  from  multiplying 


1*K 


- (1-2  € ’31+  £ 


, 9 ,-8 T.  3 ✓ — 2sT  1 , -3sT 

t - it  £ 


(68) 
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Figure  U3  shows  the  output  sanples  arising  from  an  input  which  io  a sampled 
step*  Again  the  fact  that  the  system  must  have  overshoots  to  equal 


tiaja 


rig.  k)  itesoonse  of  proposed  system  to  a step  input, 

undershoots  in  its  step  raenonse  to  have  a zero  velocity  error  coefficient 
and  the  fact  that  the  transient  is  short  requires  high  overshoot.  Assuming 
that  tii ’ >7  transient  is  acceptable,  let  us  evniiinte  K eventually 


/ * 

it  V S } • 

c 


From  F.q.  67,  K can  be  calculated  to  be 


_ -sT  . -2sT  , - 

9 r 3 / ^ i / 

I ^ -?C  • Ev 


JsT 


(i-  €“-*)"  U-6"°7U) 


'-yO- 


iK^s; 


i— 


c • 
S \ 


s+i) 


9 =* T 3 -2sT  . a _-3sT 

HC  " ? t " h ^ 

(1-  <TSV  (1-  €sl  A) 


(70) 


Given 


1C  (s)  -7- 


C -"(  3*D  } 


i ..  one  is  /Lot  forced  to  go  lr.iiquely 


to  K (s)  - out  can  go  to  a function  K having  only  Doles  within 

c s-'e  s*D  c-s* 

the  central  strip  - *^V’2  ^ cu  2.  To  this  F_c  one  can  add  any 

"zero-sampie*  K}  that  is,  one  whose  impulse  response  goes  through  zero 
at  all  sampling  points. 


Lve  lie  (s)  -n—~ I 

! c s1 2(  s+l)  I 


into  Dart,ial  fractions  so  as  to 


be  able  to  identify  terms- 


f 9 3,“?7  _ . 

MI  * * t + e c / 

(l-  C"Si)3  (l-  <fsT  A) 


1 „ ->8TV  y .-«T  . raT 

ft  > =f 

«*  ♦ — a- 


_ rn 

-Si 


7.  ,-sT\* 


u-t 


(i-  fV 


1 

27 


1 

27 


1-  c "3  r i-f'OAA 


(71) 
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The  £ (s)  --r— having  poles  within  the  etrir  -•£)■/■>  s 

3 s"(  rsvj.) 

be  identified  from  Eq.  71  by  methods  described  in  Appendix  B. 


' J' -/ 2 can 


It  I£*  2Z_.  21 

C • S • ^ Cm  _ 

ess  s+1.38? 


(72) 


i:6ls' *3.U?«*1.95 
8^(8*1.387) 


K alone  will  not  lead  to  a realizable  K (s)— because  it  has 

c*s*  c s2(r^i) 

only  two  more  poles  than  zeros.  Add  to  K a "ssro-acunyle"  K such  that 

Ci?c 

the  combination  has  3 more  poles  than  zeros.  One  such  zero-sample  K is 
given  by  Eq.  73. 


JO 

-uiL 


z.a,  s*3*  Jsr 


.UJL 


s*  3-jsr 


?«P 

2 , . . - . ' 


-2tfL 


a "*68+18.3? 


(73) 


Select  the  value  of  8 »uch  that  K * K has  3 more  poles 

c#s*  Z • 8 • 


than  zeros*  For  this  to  h-  the  case  2ffS  * 


Kc(s) 


S ( 8*1) 


1.6lB“*3.<j99*1.8*>  . -1 .61 

- — ♦ ~T, 

e-’f s+1-387)  s'*  6s*l3.89 


1C  .91$s3*5'6 . 8l?2*77.1a*3U  .93 

_ — — a 

9 (s*l .36? ) (s- ^63+18 .69) 


VB/ 


10 .9l5a3*$6 . 815>?7  .ls+34 .98 

1.3878(0^*68*18.89) 


(75) 


Hie  K^vs,1  is  certainly  physically  realizable.  The  only  trouble  with  it 


i r t.b  at 


.is  probably  more  complicated  than  it  really  neeas  to  be 
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Sucii  complication  should  have  been  anticipated  because  we  did  not  concern 

ourselves  with  realization  problems  until  the  end  of  tie  design  procedure. 

The  design  procedure  .just  discussed  selected  a desirable  — ■ 

1*K 

function  which  vas  realizable  arid  proceeded  to  derive  the  transfer  function 

of  the  compensating  filter.  Though  tue  procedure  w«s  workable  it  gave 

too  much  attention  to  system  response  and  too  little  attention  to  the  form. 

of  the  compensation.  An  alternative  procedure  would  be  to  find  what  functions 

can  be  added  to  the  fixed  S{s)  function,  so  as  to  leave  the  excess  of 

*1 

poles  over  zeros  at  3 and  still  provide  a more  desirable  — -■■■■-  . Actually 

* 1+K 

by  adding  to  K(a)  one  would  add  to  K and  it  would  be  easier  to  evaluate 
the  new  lAi£  than  its  reciprocal.  For  this  cut— and-try  design  no  neat 
schemes  have  been  worked  out  for  the  general  case  and  considerable  computation 
and  plotting  is  required. 

Multiple-loop  systems  can  be  designed  much  as  they  are  in 
conventional  servos  but  the  procedures  will  not  be  illustrated  in  this 
report ; TVi*  procedure  which  has  been  used  is  to  reduce  the  compensation 
problem  to  an  equivalent  single-loop  problem  and  to  derive  conditions  upon 
the  system  so  that  the  compensation  is  realizable.  Overall  response 
characteristics  are  chosen  and  the  compensation  network  is  derived. 

3 Summary  of  Design  Procedures 

The  first  problem  in  designing  sampled-data  servomechanisms 
is  the  problem  of  selecting  an  output  member  and  a basic  sampling  rate 
from  the  general  perfoimance  requirements  of  the  system,  Tbe  output 
member  is  selected  or.  the  same  basis  as  the  output  member  would  be  selected 
for  a continuous  servo.  The  sampling  rate  is  determined  by  the  bandwidth 
required  in  the  system  loop.  After  the  fixed  parts  of  the  system  are 
selected;  the  compensation  is  chosen « To  select  the  compensation  one 
adopts  a trial  and  error  procedure  to  approximate  the  desired  response 
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characteristics  in  terms  c?  the  compensation  filters  which  are  pnysicaily 
realizable.  The  first  part  of  this  section  dealt  with  simple  approximate 
character ization  of  the  responsos  of  a sample ci-data  system.  The  second 


part  dealt  with  physical  real  ity  ■Jcn*-iibr.or;s 


the  last  part  dealt 


with  the  realizing  the  response  characteristic:'  in  terms  of  the  framework 
of  physical  realizability  limitations.  The  Is.st  procedure  must  fundamental! 
be  cut-nnd-try  tod  a good  procedure  must  be  characterized  by  short  and 
simple  trial  step3. 
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APPENDI7  A 

Procedures  for  Evaluating  K (3)  from  K(s)  and  *6^  (s)  from tC;  ( s ) 

As  has?  been  printed  out  earlier,  the  relation  of  !f*(a)  to 
K(s)  is  the  same  as  the  relation  of  (s)  to  ^ ( - ) though  'A(s)*  and 
Kfs)  are  transfer  functions  and  (g)  and  ^(5)  are  transforms  of  time 
functions.  K (s)  and  K(s)  may  be  visualized  as  being  transforms  cf  an 


ir, pulse  response.  The  problem  discussed  in  this  appendix  will  be  t-ho 

problem  of  relating  M^*(s),  the  !*apl see;  Transform,  of  an  inpulse-modUiated 

signal,  to  0 } y the  I.aplace  Transform  cf  the  inx  sampled  sxgnal.  Because 

of  the  f?r-t  that  the  relation  be  tween  •*»,  ( s ) and  ±&.  (3)  is  linear,  suDerDcsition 

xx  * 

4 

of  component  ^ ' s leads  to  superposition  cf  KK  1 3,  and  it  is  wise  to  build 

up  c tabic  wf  pairs.  Bj  using  uhis  table  one  can  go  xroiri  a given  to 

* , 

r8^  by  resolving  one  **<*,  x«lo  oompeneno  ■ a which  are  avaxxanie  in  the 

table,  associating  each  with  ita*^,  and  finally  superposing  *9^* s. 

....  . * . 

«eeaieB*  to  say,  the  table  forms  a two-way  street  and  the  same 

4 

sort  of  procedure  can  be  used  for  associating  a^  with  a possible  *9^ . 

The  inverse  process  will  be  discussed  in  Appendix  B, 

Since  most  >9^  functions  are  rational  Lr  s.  it  is  easy  to  resolve 
W*  into  a partial  fraction  expansion.  Accordingly-  the  table  should  include 

where  0/  is  a Comdex  constant.  Table  A1 


9.  functions  of  the  form  — - — — 
i s<-o< 

below  gives  a list  of  useful  transform  pairs  between  continuous  and 


impulse -modulated  time  functions.  Because  of  X»ilG  fact  that 


corresponds  to  a function  having  a discontinuity  at  t-0  which  is  a sampling 
point,  there  is  uncertainty  about  the associated  with  it.  The  next 
few  paragraphs  will  explore  and  remove  the  confusion  arising  from  this 
uncertainty. . 
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• *1 
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Pairs  of  Transforms  of  Continuous  Tims  Functions  and 

| Their  Impulse-Mcdualted  Equivalents 

i 


Rote:  All  continuous  time  functione  having  transforms  rational  in  s 

and  having  impulse  modulated  equivalents  can  to  related  to  their 
impulse  modulator  equivalents  by  use  of  the  table  sr.d  pe.rt.iai 
fraction  expansions.  The  same  is  true  for  the  process  ox  going 
from  the  impulse -modulated  funotioi  ; to  the  equivalent  time  function* 


i 


i 


Report  R-222 


*$- 


According  to  the  view  6f  the  impulse  modulator  os  a sideband  generator, 

Stag- 


^ >|  \ 


wnere 


«*<c 


f 4.0' 


?*• 


T «*  4 v%+ » 


OS  uWS«ii  C<3?if“>JL05  i 


X # 

If  #©t(s)  *» . the  value  oT&{  (s)  can  be  evaluated  as 

I 3*  er 


"I  w • $ > 


n-- < 


s-<  + jn  A. 


(41) 


Equation  (Al)  can  be  verified  by  a summation  process.  This 
result  agrees  with  the  one  given  in  Eq.  5 of  Section  .1  tin  the  text. 

The  transform  — — r <v>»re  spends  to  s string  of  impulses  in  the  rime  domain 
whose  successive  areas  are  v"‘*  with  a first  impulse  whose  area  is  ^ 

because  only  half  ox  the  first  sample  is  caught  by  an  ideal  impulse 
modulator.  The  value  of  *©l  can  be  obtained  by  adding  up  the  transforms 


of  one  individual  impulses. 


*r-sT  ^ ~2*T  * »2aT 

€ * r s r 


* 


1 

~ £ )T 

(A2) 


1 

2 


The  equivalence  cf  Eqs.  (Al)  and  (A2)  indicates  agreement  between 

the  "complementary  signal"  approach  and  the  "shifted  impulse"  approach. 

Which  of  the  two  one  should  use  depends  upon  which  summation  j.3  easier 

to  form.  The  fact  that  the  time  function  corresponding  to — - — has 

Sv  «< 

a discontinuity  at  the  origin  which  is  a sampling  point  and  one  might 
not  wish  to  assume  that-  physical  equipment  would  catch  only  half  of  the 
first  sample.  Accordingly  one  can  advance  the  wave  by  shifting  it  to 
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the  left  by  a small  amount  r(f  . The  first  sample  caught  then  would  become 

Tl 

1 instead  of  and  psrfc-'-»-ir»g  the  approach  of  Eq.  (A2)  >?ould  yl old 


(a) 


1 . 

rzn — rvm 

/•.  -VO’  Of.  t J- 


(A  i) 


A check  with  the  approach  of  Eq.  (Al)  e«n  be  obtained  by  performing  the 
summation  of  all  complementary  signal  components. 


***(«> 


w>0 

11m  1 \ 

* * > 


(• 


S(s+jn-fi) 


+ .1n.fi 


7"-=r8+3STr 


i >i<  ) 

• \ '*H/ 


If  none  of  the  first  sample  were  caught,  then  the  summation  of 


shifted  impulses  would  start  with  the  one  which  occuri cd  at  t-T,  the  transform 


of  which  is  (~ 


-oT 


- -J  ? 


x ^ Tiie  approach  of  Eq.  (k2)  would  yield 


£-(s+°<)T 

' / dsWIT 

L ‘ 


TT — rtir 


(A5) 


To  insure  that  the  sample  «t  t-0  is  net  caught,  the  time  function 
could  be  delayed  by  £ seconds  where  £ is  a very  small  number.  The 
approach  of  Eq.  (Al)  can  be  checked  by  again  performing  a summation  on 
all  complementary  signal  components. 


<M  - . I L_ 


f c £ i ^ ^ 


o(  .jn  JX 


: — :r?oir:/Tf  ~ J 

t w *A  • 

(A6) 


Summarising  the  results  of  Eq3.  (Al)  through  (A6),  the  relation 

1 + 

between  — -- — y and  the  corre  ■<..  -rdijig  ?0t  depends  upon  what  assumption  is 

a*o(,  - 

made  about  the  size  of  the  sample  taken  at  the  discontinuity.  The  choice 
made  in  Table  Al  was  made  arbitrarily  and  only  because  of  the  simplicity 
Jn  the  resulting  formulae.  In  practical’,  cases,  few  sampled  time  functions 
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really  have  discontinuities  at  the  origin  and  accordingly  the  sum  of  the 
residues  is  zero.  This  fact  means  that  which  form  is  chosen  is  inconsequential , 
As  .an  example,  assume  F(s)  is  given  by  Zq.  (A7). 

*1 


m — «=»«*—  A 

S+o>1  s+  o<2  s+a<f3 


where 


,al  + a2  + a3  ' °* 

/lake  the  identification  of  E a.  (A8). 


v s*  < 


. ; _ 


1 


1-  e 


♦ p. 


(A?) 


(A8) 


f 1.0  , — ~ , — 

• i.  ^ -u*  srpr 


a.. 

.1 


+ |jrj . V 


l-  5 

* * *y 


(A  9) 


F (s)  does  not  involve  JJ  since  the  sum  of  its  coefficients  add 
to  “ere-  Thin  fact  viii  always  be  the  case  in  practical  situations  because 
no  physical  quantities  really  are  represented  by  functions  having 


disc 
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APPENDIX  B 

j£. 

PpnroHn res  f or  Evaluating  KCs)  from  K (s)  audit's)  from  (a) 

The  problem  in  going  from  to  iV^(s)  has  really  one  more 

, •Q-/  . 

complication  than  tlie  problem  of  going  from^ta)  to  (s)  because  the 
continuous  function  cannot  be  uniquely  defined  from  its  samples.  The 
procedure  arbitrarily  adopted  here!  is  to  select  the  iw^(b)  having  its 
poles  in  the  strip  around  the  origin  in  the  s -piano . This  A (s)  is 
referred  to  as  the  "central  strip" >0^.  One  can  add  to  this >©<  c g 
any  transform  corresponding  to  a time  function  whose  samples  all  are  zero 
Such  a function  need  not  be  analytically  defined  but  in  practical  cases, 
particularly  in  connection  with  relating  sampled  and  continuous -signal 
transfer  functions,  it  car.  be  most  conveniently  used  if  the  having 
zero  samples  is  rational  ir.  s.  Two  kinds  of  "zero  sample"  ' s have 

been  used.  One  kind  is  made  up  of  functions  whose  poles  lie  on  the 
lines  s » ♦jk.fi /2  and  whose  residues  in  these  poles  are  imaginary. 

Each  of  these  functions  will  be  aasoeiatod  with  tirsc  functions  represente 
by  rotating  vectors  having  zero  res!  Dart  at  the  sampling  instants. 

The  other  kind  is  made  up  of  functions  having  poles  separated  by 
in  the  s-plane  and  having  residues  whose  sums  are  zero.  As  illustrated 
in  the  exanples  in  the  report,  it  is  often  possible  to  use  the  freedom 
from  non-uniqueness  in  going  fromi^  (s)  to  t©^(s)  to  practical  advantage 
but  full  advantages  of  the  use  of  this  freedom  have  not  bets-,-  colored* 


